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Abstract. We present pointwise gradient bounds for solutions to p-Laplacean 
type non-homogeneous equations employing non-linear Wolff typc potcntials, 
and then prove similar bounds, via suitablo calorie potentials, for solutions to 
parabolic equations. 
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1. Introduction and main results 

The aim of this paper is to give a somehow unexpected but nevertheless nat- 
urai maximal order - and parabolic - version of a by now classical result due 
to Kilpelàinen & Maly [35], later extended, by mcan of a different approach, by 
Trudinger & Wang [55] . In [351 |5S] the authors obtained a neat pointwise bound 
- see (|1.14p below - for solutions to non-homogeneous quasi-linear equations of p- 
Laplacean type by mean of certain naturai Wolff type non-linear potentials of the 
right-hand side measure datum; this is commonly considered as a basic result in 
the theory of quasi-linear equation. In this paper we upgrade such a result by show- 
ing that similar pointwise bounds also hold for the gradient of solutions, scc (jl.6p 
below. The resulting estimate finally allows to draw a complete and unifled picture 
concerning gradient estimates for degenerate problems: it permits to recast clas- 
sical gradient estimates for the non-homogeneous p-Laplacean equation obtained 
by Iwaniec [29] and DiBcnedctto & Manfredi [18], as well as the pointwise L°° 
estimates obtainable whcn the right hand side is integrable enough [45] [46] , up to 
the classical gradient bounds of Boccardo & Gallouèt [TU] E] , Lindqvist [U and 
Dolzmann & Hungerbiihler & Miiller [19| for measure data problems; for more on 
the p-Laplacean see also the notes of Lindqvist [33] • More significantly, the results 
here allow for a characterization of Lipschitz continuity of solutions in terms of 
Riesz potentials which is analogous to the one available for the standard Poisson 
equation; see (|1.8p and p.l6p below. Furthermore, the estimates extend to the case 
of coefficient dependent equations, allowing to get regularity criteria with respect 
to the regularity of the coefficients which are exactly those ones suggcsted by the 
analysis of the fundamental solutions of linear equations [55] . Our elliptic results 
find a naturai but non-trivial analogue in the parabolic case, as we are also going to 
show here for a significant class of parabolic problems. In this case we shall employ 
a naturai family of "calorie" type Riesz potentials to give a pointwise bound for 
the spatial derivative of solutions; in turn, the resulting inequality implies relevant 
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sharp borderline integrai estimates which, known in the elliptic case, are otherwise 
unreachable via the known techniques, in the paraboUc one. In particular, this 
finally allows to make the missing link between elliptic and parabolic borderline 
estimates such as those in Lorcntz spaces. Moreover, as in the elliptic case, we 
find a sharp characterization of the locai boundedness of the gradient via Calorie 
potentials] see (jl.27p below. It is worth mentioning that the Lipschitz continuity 
estimates and the techniques of this paper are the starting point for furthcr de- 
velopments, eventually leading to C^-regularity assertions for solutions [53]. For 
notation and more dctails wc recommend the reader to look at Section 2 below. 

1.1. Degenerate Elliptic estimates. In this section the growth exponcnt p will 
be a number such that p > 2, whilc for the rcst of the paper O will denote a bounded 
open and Lipschitz domain of R", with n > 2. We shall consider general non-linear, 
possibly degenerate elliptic equations with p-growth of the type 

(1.1) — div a{x, Du) = /i . 

whencver is a Radon measure defined on fi with finite total mass; eventually 
letting /i(IR" \ ri) = 0, without loss of gencrality wc may assume that /i is defined 
on the whole M". The continuous vector field a: fi x R" — > R" is assumed to be 
C^-regular in the gradient variable z, with a^^-) being Carathéodory regular, and 
satisfying the foUowing growth, ellipticity and eontinuity assumptions: 

r \a{x,z)\ + \a,{x,z)\i\z\^ + s')i <L{\z\^ + s^)^ 

(1-2) I i.-i(|z|2 + s2)^|A|2<(a,(x,z)A,A) 

[ \a{x,z) - a{xo,z)\ < Liuj{\x - xo\){\z\'^ + s^)^- 

whenever x,xo G il and z,\ G R", where < u < L and s > 0, Li > 1 are fixed 
parametcrs. In p.2p the function uj : [0, oo) [0, oo) is a modulus of continuity i.e., 
a non-decreasing concave function such that wiO) = = luHgiQLo^g) and lu{-) < 1. 
On such a function we impose a naturai decay property, which is essentially optimal 
for the result we are going to have, and prescribes a Dini-eontinuous dependenee of 
the partial map x i— > a{x, z)/{\z\ + s)p~^: 

(1.3) rHg)]-p^:^d{R)<oo. 

Jo Q 

The prototype of (jl.ip is - choosing s = and omitting the x-dependence - clearly 
given by the p-Laplacean equation 

(1.4) ^àÌY{\Du\P-^Du)^ ^i. 

Our estimate involves the classical non-linear Wolff potential defined by 

The first result, which we naturally state as an a priori estimate, is 

Theorem 1.1 (Non-linear potential gradient bound). Let u G C^{n) be a weak 
solution to (|l.ip with ^ G L^{ri) under the assumptions (jl.2p . Then there exists a 
Constant c = c{n,p, v, L) > 0, and a positive radius R = R{n,p, v, L, Li, U!{-)), such 
that the pointwise estimate 

(1.6) \Du{xo)\<cJ {\Du\ + s) dx + cWi {xo,2R) 

Jb{xo.R) p'^ 

holds whenever B{xq,2R) C iì and R < R. Moreover, when the vector field a(-) is 
independent of x - and in particular for the p-Laplacean operator (jl.4p - estimate 
(|1.6p holds without any restriction on R. 
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The potential appearing in (|1.6p is the naturai one since its shape respects 

the scaling properties of the equation with respect to the estimate in question - see 
Section rOl and compare with the linear estimate in (|1.12p : moreover, the potential 
is in a sense optimal, as in an estimate like (|1.6|) it implies the recovering of 

the sharp integrability results known for general weak solutions to (|1.1[) . Moreover, 
as a matter of fact when the right hand side is a properly singular measure estimate 
p.6p reverses; for ali such aspects we refcr to Rcmark l6.2l below. An approximation 
procedure cventually allows to remove the extra regularity assumptions u € C^(f2) 
and ^ e L^(ri) up to consider the most general case of solutions to measure data 
problcms as 



(1.7) 



— div a(a;, Du) = /.i in 

u — {) on do. . 



Our result indeed holds for general Solutions Obtained by Limit of Approximations 
(SOLA) to (|1.7p . a class of very weak solutions - solutions not necessarily lying 
in the naturai energy space W^'P{^) - which is introduced to get existence and 
uniqucness in several cases. Moreover, SOLA coincide with usuai weak solutions 
for rcgular data /i S W~^'P (lì); we refer to ScctionObclow for more details. see also 
Thcorcms l5.Hl5^ bclow. When extended to general weak solutions estimate (|1.6p 
tells US the remarkable fact that the boundedness of Du at a point xq is independent 
of the solution u, and of the vector field a(-), but only depends on the behavior of 
1^1 in a neighborhood of xq. 

Corollary 1.1 (C°'^ -regularity criterium). Let u e W^-P^\fl) be a SOLA to the 

problem - which is unique in the case ^ G L^(iì) - under the assumptions 

(HHl) and do]). Then 

(1.8) {■,R) e L'^in) for some R>0=> Due L^^{n,m:"). 

Furthermore, there exists a Constant c, depending only on n, i/, L, Li, d{-), such that 
the following estimate holds whenever B2R C CI: 



(1.9) <c/ {\Du\ + s)dx + c {■,R) 



Remark 1.1. CoroUarv 11.11 allows for a naturai Lipschitz continui ty criterium 
with respect to the regularity of coefìicients (|1.3|) and moreover to obtain for the 
p-Laplacean, cssentially the samc criterium of Lipschitz continuity availablc for the 
Laplacean, that is the one via Ricsz potentials; sce (jl.l6p below and also (|6.4p for 
a corollary. We remark that finding conditions on jjL implying the boundedness of 
Du was a major open issue in the theory of p- Laplacean equation and p.Sp seems 
to provide a first satisfying answer to the problem. Indeed Cianchi |13| has shown 
that already in the plain case of the standard Poisson equation 

(1.10) -Au = ^i 

the gradient might be unbounded as soon as 2 ^ The estimates above also 
emphasize the role played by Dini-continuity of the coefficients in order to get the 
boundedness of the gradient, which is basically the same one observed when dealing 
with linear elliptic equations with variable Dini-continuous coefficients, and proving 
pointwise bounds on the derivatives of the Green's function [25l Section 3] . As a 
matter of fact a feature of the approach adopted here is that estimate (|1.6p allows 
to derive conditions for gradient boundedness, which are borderline simultaneously 
w.r.t. the right hand side - i.e. (jl.Sp - and w.r.t. coefficients - i.e. (|1.3p . already in the 
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linear case. Indeed, a recent striking example of Jin & Mazya & Van Schaftingen 
[50] shows that there exist distributional solutions u to linear equations of the type 

div {A{x)Du) = 

with a continuous (not Dini-continuous) and elliptic matrix A{-), such that u G W^'P 
for every p < oo as prcdicted by the standard rcgularity thcory. but such that 
u ^ W'^'°°; in this case one evcn has that Du ^ BMO. 

1.2. Zero-order and related estimates. The linkage between estimate (|1.6p and 
the estimates of [35l [56] is clear. The hnear case a{x, z) — z obviously Icads to the 
standard Poisson equation (|1.10p - here for simplicity considered in the whole M" - 
for which, due to the use of classical representation formulas, it is possible to get 
pointwise bounds for solutions via the use of Riezs potentials 

(1.11) /,(^)(,):^£_M|L_, ^g(o,n] 
such as 

(1.12) \u{x)\ < cl2M){x) and \Du{x)\ < ch{\^i\){x) , 

the first being actuaUy valid for ti > 3. We recali that the equivalent, localized 
version of the Riesz potential Ip{fJ.) is given by the linear potential 

(1.13) r^ixo, R) ^ Wl^ixo, R) - £ ^^^^^qr^ j < , 

the last line being valid for non-negative measures (note that we use two different 
- but similar - notations for the Riesz potential in (|l.lip and its "polar version" in 
(|1.13p ì. The naturai non-linear version of (|1.12p . derived in |351I56| for non-negative 
measures, when p < n, is treated via Wolff potentials: 



(£ 



(1.14) \u{xo)\ < c [ -I- {\u\+s)^dx\ +cW^p(xo,2i?) 7>p-l. 

lB(xo.R) J 

Estimate (|1.6p upgrades (|1.14p to the gradient/maximal level, obviously replac- 

ing with , and represents the p-Laplacean analog of the second es- 

p 'P 

timate in the left hand side of (|1.12p . Indeed, for p = 2 by (|1.13p we have 

W'^ (xo,i?) < /i(|^|)(xo). Let US mention that the technique developed for esti- 

2 

mate (|1.6p also yields an alternative proof of estimate (|1.14p which now holds for 

general signed measures; see Remark 14.11 below. We also refer to [49] for gradient 

potential estimates when p = 2, and to the important work of Labutin [35] for 

relevant Wolff type potential estimates related to fuUy non-linear Hessian type op- 

erators. For further interesting relations between degenerate quasilinear equations 

and potentials we refer to the recent interesting paper of Lindqvist & Manfredi [H] . 

Another consequence of estimate (|1.6p and of the classical bound 

p'P Jo V Q J Q 

(1.15) < c/Jk(|H)l^j(xo)=:cV,.^(|H)(xo) 
is the estimate 

(1.16) \Du{xo)\<c~f {\Du\ + s)dx + cVij,{\fx\){xo) 

which holds whenever B{xo,2R) C il satisfies the smallness condition imposed 
in Thcorcm 11.11 Here we remind the reader that we have previously extended n 
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to the whole space R". The non-hnear potential Vi (/x)(a;o) - often called the 
Havin-Maz'ja potential of /i - is a classical objcct in non-linear potential theory, 
and together with the bound (|1.15p comes from the fundamental and pioncering 
work of Adams & Meyers and Havin & Maz'ja; see also [HEll^ZI- Estimate (jl.l6p 
allows to derive ali types of locai estimates starting by the properties of the Ricsz 
potential; see Section [6] below. 

1.3. Parabolic estimates. Our aim here is not only to give a parabolic version 
of the elliptic estimate (|1.6p . but also to give a zero order estimate, that is the 
parabolic analog of the zero order elliptic estimate [35j . which at this point will 
essentially foUow as a corollary of the proof of the gradient estimate. We consider 
quasi-linear parabolic equations of the type 

(1.17) Ut — dìv a{x,t, Du) = ^ , 

in a cylindrical domain J7t := £1 x (— T, 0), where as in the previous sections f2 C 
M", n > 2 and T > 0. The vector-field a: x M" ^ R" is assumed to be 
Carathèodory regular together with az{-), and indeed being C^-regular with respect 
to the gradient variable z £ M", and satisfying the foUowing standard growth, 
ellipticity/parabolicity and continuity conditions; 

r \a{x,t,z)\ + \a^{x,t,z)\{\z\ + s) < L(\z\ + s) 

(1.18) < iy\X\^ < {a^{x,t,z)X,X) 

{ \a{x, t, z) — a{xo,t, z)\ < Lilu{\x — xo\){\z\ + s) 

for every choice of x,xo e fi, z,X e M" and t E (— T, 0); here the fmrction 
uj: [0,oo) [0, 1] is as in ([I2])3 for p ~ 2. Note that anyway we are assuming 
no continuity on the map t i—t a(-,i, •), which is considered to be a priori only 
measurable. In other words we are considering the analog of assumptions (jl.2p for 
p = 2; the reason we are adopting this restriction is that when dealing with the 
evolutionary p-Laplacean operator estimates take the usuai forni only when using 
so called "intrinsic cylinders" , according the by now classical parabolic p-Laplacean 
theory developed by DiBenedetto [Ij . These are - unless p = 2 when they reduce 
to the standard parabolic ones - cylinders whose size locally depends on the size of 
the Solutions itself, therefore a formulation of the estimates via non-linear poten- 
tials - whose definition is built essentially using a standard family of balls and it 
is therefore "universal" - is not clear and will be the object of future investigation. 
We refer to [2] for global gradient estimates. 

In order to state our results we need some additional terminology. Let us recali 
that given points {x,t), {xo,to) £ R"+^ the standard parabolic metric is defìned by 

(1.19) dpi,r{{x,t),{xo,to)) := max{|a; - xo\, ^y\t - to\} « ^y\x - xqP + \t - to\ 

and the related metric balls with radius R with respect to this metric are given by 
the cylinders of the type B{xo, R) x {to~ R'^ ,to + R^). The "calorie" Riesz potential 
- compare with elliptic one defined in p.lip and see Remark 11.21 below - is now 
built starting from 

(1.20) I,{^.){{x,t))■.^ I frfi^N-p 0<P<N:=n + 2 

whenevcr {x,t) G R"+^. In order to be used in estimates for parabolic equations, 
it is convenient to introduce its locai version via the usuai backward parabolic 
cylinders - with "vertex" at (xo,to) - that is 

(1.21) Q{xo, to; R) B{xo, R) x (io - i?', , 
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and is now given by - with N :~ n + 2 

(1.22) I,1(.o,^o;i?):= rM%|£i^^!£ where/3e(0,iV]. 
The main rcsult in the parabohc case is 

Theorem 1.2 (Calorie potcntial gradicnt bound). Under the assumptions (|1.18p 
and (fO)) . let u <E C°(-T, 0; ^^(O)) be a weak solution to PT7)) with fi G L'^{Q.t), 
and such that Du G G'^{^t)- Then there exists a Constant c = c{n,i',L) and a 
radius R = iì(n, v, L, Li, a;(-)) such that the following estimate: 

(1.23) \Du{xQ,tQ)\<ci {\Du\ + s)dxdt + cmxo,to;2R) 

JQ(xo,to;R) 

holds whenever Q{xo,to',2R) C fi is a backward parabolic cylinder with vertex at 
{xo,to) and such that R < R. 

As in the ehiptic case when deahng with parabohc equations with no space 
variable dependence of the type 

(1.24) Ut - dÌY a(t, Du) = fi 
we can improve the result of Theorem 11.21 as foUows: 

Theorem 1.3 (Componentwise calorie bound). Under the assumptions (jl.lSp . let 

u G C''(— T, 0; i^(ri)) be a weak solution to <\\.2A\ with fi G ì^(ÌÌt), and such that 
Du G C^{VIt). Then there exists a Constant c = c(n,v,L) such that the following 
estimate: 

(1.25) \D^u{xo,to)\ < cj- {\D^u\+ s)dxdt + cI'^{xo,to]2R) 

jQ(xo,to;R) 

holds whenever Q{xQ,tQ;2R) C fi is a backward parabolic cylinder with vertex at 
{xo,to), and ^ G {1, . . . 

Estimates (|1.23p - (|1.25p also hold for general weak and very weak solutions, and 
in particular for solutions to measure data problems as 



(1.26) 



Ut — diva(x', t, Du) ~ fi in ^It 

U = On dpur^T : 



where fi is a general Radon measure on CIt with finite mass on CIt, that we shall 
consider to be defined in the whole R"+^. See next section for the definition of 
solutions. In the spirit of the elliptic result (|1.9p we have the following implication, 
whieh provides a boundedness criteria for the spatial gradient, under the Dini- 
continuity assumption for the spatial coefficients ([TT 



(1.27) Ii(-;^) e L'^inr), for some R>0=^Due L^^{nT,W) . 

We conclude with the zero order potcntial estimate, which applies to general equa- 
tions of the type (jl.lZp when considered with a measurable dependence upon the 
coefficients {x,t). The relevant hypotheses in this case are the following standard 
growth and monotonicity properties: 



(1.28) 



\a{x,t,z)\<L{\z\ + s) 
v\z2 - Zip < {a{x, t, Z2) - a{x,t, zi),Z2 - zi) 



whieh are assumed to hold whenever {x,t) G flr and z,zi,Z2 G M". In particular, 
since the pointwise bound will be derived on u, rather than on Du, we do not need 
differentiability assumptions for a(-) with respect to the spatial gradient variable z. 
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Theorem 1.4 (Zero order estimate). Under the assumptions (|1.28p . lei 

u e L^{-T, 0; W^^'^in)) n C"(17t) 



be a weak solution to (jl.l7p wit/i /i S ^^(Oj'). T/ien i/iere exists a Constant c, 
depending only on n,v, L, Li, such that the foUowing inequality holds whenever 
Q{xQ,tQ;2R) Q is a backward parabolic cylinder with vertex at (xo^to): 



Remark 1.2. The calorie Riesz potential defincd in p.20p is difFerent from the 
parabolic Riesz potential considered for instance in [3] via convolution of /z with 
the heat kernel 



On the other hand both the kernels considered exhibit the sanie scaling with respect 
to the parabolic dilation R — s- (Rx, R'^t), which is in turn the relcvant property to 
determine the regularization properties of the related convolution operator. There- 
fore, as we shall see from Section lOl below. we have that the calorie Riesz potential 
p.20p is perfectly sufficient to infer the sharp regularity/integrability properties 
of Solutions. In particular, from (|1.23p we shall be able to derive a borderline 
Marcinkicvicz estimate which seems to be difhcult to derive using for instance the 
truncation methods from |n|. Anyway. for estimatcs regarding (|1.30p see [57] . 

1.4. Pian of the paper. After establishing in Section 2 some notation, Section[3] 
is dedicated to the proof of Theorem 11.11 This will require a carcful combination 
of regularity estimates for p-harmonic functions, together with suitable comparison 
estimates which make the density of the Wolff potential appear; ali such ingredients 
will be finally combined in a suitable iteration scheme. The same scheme will be 
foUowed in Section U] for the parabolic case, where additional difficulties come into 
the play; in particular, the proof of the comparison estimates necessitates further 
delicate technicalities, whilc some precise estimates from Nash-Moser's theory will 
be needed. In Section 5 we show how the cverywhcre pointwisc estimates derived 
for a priori regular solutions, actually extend to ali kinds of general weak solutions, 
including solutions to mcasure data problems. which are the most general ones. 
Finally, in Section 6 we demonstrate how the pointwisc estimatcs derived here 
allows to recast essentially ali the main gradient estimates for non-homogeneous 
problems as described at the beginning of this Introduction, and in particular those 
in finer scales of spaces such as Lorentz or Orlicz spaces. We re-emphasize here 
that, when considering parabolic problems, the techniques presented here are the 
only one available for getting such estimates. The techniques presented in this 
paper are general enough to be applied in different contexts; an example is |12j . 
where problems with non-standard growth of p(x)-type - see for instance [T] - are 
considered. 

Some of the rcsults here have been announccd in the nota lincea |22j . 

Acknowledgments. This research is supportcd by the ERG grant 207573 "Vec- 
torial Problems" , and by MIUR via the national project "Calcolo delle Variazioni" . 
We also thank Verena Bògelein who carefuUy read a draft of the paper. 



(1.29) 




(1.30) 





2. NOTATIONS 



In this paper we foUow the usuai convention of dcnoting by c a general Constant 
larger (or equal) than one, possibly varying from line to line; special occurrences 
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will be denoted by ci etc; relevant dependence on parameters will be emphasized 
using parentheses. We shall denote in a standard way 

B{xo,R) := {x e M" : \x - xo\ < R} 

the open ball with center xq and radius i? > 0, wliile backward parabolic cylinders 
have been already defined in (jl.2ip : when not important, or clear from the context, 
we shall omit denoting the center as foUows: Bji = B{xo,R), and the sanie will 
happen for parabolic cylinders concerning the vertex. Moreover, when more than 
one ball - resp. cylindcr - will come into the play, they will always share the same 
center - rcsp. vertex - unless otherwise stated. We shall also denoto B = Bi = 
i3(0,l), and Q = Qi — Bi x (—1,0). We recali that, given a cylindrical domain 
of the type C ~ Ax (ti, ^2) - and in particular a parabohc cy Under - its parabolic 
boundary 9parC is given by dC\ {Ax {^2})- With A being a measurable subset with 
positive measure, and with g : A — > R'^ being a measurable map, we shall denote 

g{x)dx -i- / g{x)dx 
A 1^1 Ja 

its average. In this paper, ali the measures considered will be Radon measure with 
finite total mass, for convenience assumed to be defined on the whole space: M" 
for elliptic problems, W^'^^ for parabolic ones; with A C K" or A C R""*"^, we shall 
denote by M. {A) the space of ali Radon measures with finite total mass defined on 
A. When a measure fjL actually turns out to be an L^-function, we shall use the 
standard notation, whcncvcr A is a measurable set on which /i is defined 

ImI(A) := / \tx{x)\dx. 
Ja 

When dealing with general measure data problems as in (|l.ip and p.l7p . we shall 
consider so called very weak solutions, i.e. those solutions who are not ncccssarily 
lying in the naturai energy spaces associated to such equations growth assump- 
tions ([TlDi and ([TT^ i being given - spccifically W^^P{n) and ^^(-T, 0; W'^'^{VÌ)), 
respectively - but that are nevertheless integrable enough to allow for a usuai distri- 
butional formulation. Such solutions are not necessarily unique, but suitable rein- 
forced notions of solutions can be considered as well, in order to achieve uniqueness 
in several cases; see Section 5 below. In the elliptic case a (very) weak solution to 
(|l.ip is a function u £ W^'P^^{fl) such that the distributional relation 

{a{x, Du), Dip) dx = / pdfx 
n Jn 

holds whenever p G C^{iì) has compact support in fi. In the parabolic case a 
(very) weak solution is a function u E L^{—T, 0; W^'^{n)) such that 

(2.1) — / uiptdxdt+ / {a{x,t, Du), D(p) dx dt = / (pdfi 

holds whenever ip S C^{p.T) has compact support in VIt- When referring to 
Cauchy-Dirichlet problems of the type p.26p , whilc the lateral boundary condition 
can be formulated by prescribing the inclusion u e L^(— T, 0; Wg^'^(ri)), the initial 
boundary onc u{x, —T) = is undcrstood in the L^-sense, that is prescribing that 



1 



-T+h 



lini— / / \u(x,t)\dxdt ~ Q. 

hiQ h J_T Jn 
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In the parabolic setting a convenient "slicewise" reformulation of (j2.ip is given by 
mean of so-called Steklov averages; in fact, for h > and t E [—T,0) we define 

r 1 /•*+'' - ~ 

(2.2) uhix,t):= \ h ], "^"^'^^"^^ t + h<0 

[ ìf t + h>0. 

With such a notation wc have - see [T71 Chaptcr 2] - that whcn fi £ L^{ì1t) the 
slicewise equality 

(2.3) J (dtUh^p + {[a{-,t,Du)]h,Dtp)^ dx j ipfihdx 

holds whenever G C§°{Q) has compact support in fi, and for a. e. t e (— T, 0). 
As the reader will recognize we use both the notations Wt and dtw for tlic time 
derivative of a function w. 



3. ELLIPTIC ESTIMATES 

This scction contains the proof of estimate (|1.6|1 . therefore we shall argue under 
the assumptions of Thcorcm 11.11 whilc u E C^(f7) denotes the solution identified 
in Theorem ll.il 

3.1. Basic preliminaries. Lct us recali a few basic strici monotonicity properties 
of the vcctor ficld z i-^ a{-, z) under the assumption (|1.2p 9: for an absolute Constant 
c = c{n,p, v, L) > 1 the incquality With s > bcing the number appcaring in p.2p . 
wc define 

(3.1) c~^(|zip + |22p + s^)'^|z2-ziP < {a{x,Z2)- a{x,Zi),Z2- Zi) 

holds whenever 2;i,Z2,M" and x E fi. In particular since p > 2 the previous in- 
cquality implies Indeed using ([1.2^ 9 is standard to see that the inequality 

(3.2) c~^\z2 - zi\P < {a{x, Z2) - a{x, Zi),Z2 ~ zi) . 

The next result on Reverse Hòlder type inequalities will be useful in the sequel. 
Lemma 3.1. Let g : A—t K.*^ be a integrahle map such that 

U Iffl^^'MxV" <cj \g\dx 

holds whenever B2R Q A, where A C R" is an open subset, and xa > 1; c > 0. 
Then, for every t € (0, 1] and x G (0, Xo] there exists a Constant cq = Co{n, c, t) such 
that, for every B2R ^ A it holds that 

(3.3) (l Isl^d^V <co f/ \g\'dx 
\Jbr J \Jb2r 

The assertion also holds when A C M"+^, considering backward parabolic cylinders 
of the type in (|1.2ip . instead of standard balls. 

3.2. A decay estimate below the naturai growth exponent. Here we prove 
a decay estimate for solutions to ccrtain homogcncous equations which differs from 
more standard ones in that the exponents involved are smaller than those typically 
used - the "naturai growth exponents"; this estimate could be of its own interest. 
Let US consider an energy solution v € W^-'P{A) to the homogeneous equation 

(3.4) diva{xQ,Dv)^0, 

where xq € ^ and A C fi is a sub-domain of fi, in other words we consider the 
vector field from Theorcm ll.il "frozen" at a point. 
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Lemma 3.2. Let w € Vl^^'^(A) be a solution of the linear equation of the type 

div{à{x)Dw) =0, 

where the matrix à{x) has measurable entries and satisfies the following elipticity 
and growth bounds, far every A G R"; 

c*|A|^ < {à{x)\, A) , \à{x)\ < e, , for some c* € (0, 1) and c,, > 1 . 

Then there exists constants c > 1 and G (0, 1], both depending only on n^c^^jc^, 
such that the following estimate holds whenever Bg C Bfl C A are concentric balls: 

(3.5) \w-{w)B,\dx<c(^^ \w~{w)Ba\dx. 

Proof. This rcsult is a rathcr standard consequence of DeGiorgi's theory for linear 
elliptic equations. The statement foUows directly from [321 CoroUary 1,5] or by the 
proof of Theorem 7.7], see also (7.45)], taking into account the specifìc 
equations we are deaUng with; in particular we may take x = in |24j . □ 

Theorem 3.1. Let v e W^''''{A) be a weak solution to (j3.4p under the assumptions 
p.2p . Then there exist constants (3 e (0,1] and c > 1, both depending only on 
n,p,i',L, such that the estimate 



(3.6) / \Dv~{Dv)B,\dx<c(^Y j \Dv ^ {Dv) 



B„\dx 

Br, 



holds whenever Bg C Bji C A are concentric balls. 

Proof. We shall use some iteration techniques [18l [39l [21]; in particular, we shall 
use the idcas in the important paper [39| . The proof is divided in several stcps. In 
the following B{xo,R) = B^ is a fixed ball as considered in the statement of the 
Theorem. We recali that standard regularity gives that Dv e L^^{A)] moreover 
without loss of generality we can assume that Br (s A so that in the following we 
bave that Dv € L°°{Bji). From now on, for the rest of the proof, ali the balls 
considered will bave the same center xq and r will denoto a positive radius such 
that r < R. Accordingly, we set 

E{r) := / \Dv- (Dv) B^dx, M{r) := max sup \D(^v\ . 

JBr !<?<" B^ 

Step 1: Regularization. In the following wc shall considcr the non-degenerate 
case s > 0, but we shall prove estimates which do not depend on the parameter s. 
The case s = can be deduced by approximating originai solutions by solutions to 
equations of similar type satisfying s > 0, via a completely standard approximation 
procedure; see for instance [TOl 

Step 2: The fundamental alternative. This goes as follows - see [Tni[ini[S2: With 
B2r C Bu and iiq G (0, 1), clearly one of the following three alternatives must hold: 

\{D^v < M(2r)/2} n B2r\ < /iol^srl for some ^ 

\{D^v > -Af(2r)/2} n B2r\ < tio\B2r\ for some ^ 

\{D^v <M{2r)/2}nB2r\, 

\{D^v > -A/(2r)/2} n Bari > Mol^srl for ali C ■ 

The cruciai and well-established regularity property of solutions - see for instance 
|16L I39L H5] - is that there exist universal numbers fj,Q,ri = fiQ,r]{n,p,i', L) g (0,1) - 
i.e. independent of the vector field considered a{xQ, z) and of the solution v - such 
that if one of the first two alternatives in (13.711 holds then 



(3.7) 



(3.8) \Dv\> M{2r)/4: in Br 
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otherwise, if the third possibility from (|3.7p occurs, we have that 

(3.9) A/(r) < riM{2r) . 

Ali in ali wc havc that cithcr (|3.8p or p.9p happcns. The proof now consists of 
showing that combining thcsc two alternatives (|3.6p foUows. Wc recali that under 
the present assumptions we have that v € W^io'c (^) each componcnt D^v = w 
satisfies the following difFcrcntiatcd equation: 

(3.10) div (à(a;)_Dw) = ài^j{x) := (ai) ^jixo, Dv{x)) . 

Moreovcr, standard a priori estimates for p-Laplacean type equations as those con- 
sidcrcd in (|3.4p provide a Constant Cg = Cg{n,p, L) such that 

(3.11) sup \Dv\ < Cg -t {\Dv\ + s) dx holds for every i?2r ^ Br . 

Br J B2r 

Step 3: The non- degenerate case I. In this case we assume that s > 2y/nAl{r) 
for a certain radius r < R; thcn the matrix à{x) in p.lOp satisfies on Br the bounds 

i^-^sP-^lXl"^ < {à{x)X, A) , \à{x)\ < L2P'^sP-^ for every A e M" . 

Thcrcfore the matrix à(a:;) defincd in p.lOp satisfies the assumptions of Lemma l3.2L 
that applicd te every componcnt D^v = w oi the gradient gives 



(3.12) 



/ \Dv-{Dv)B^\dx<Cd(-Y" -l \Dv-{Dv)BAdx 

J Bn ^ ^ Br 



where Cd > 1, /3o G (0, 1) depend only n,p, v, L. 

Step 4- Choice of the constante. We first take a positive Hi € N such that 

(3.13) SV^Cgf]"^-^ < l . 

This fixes Hi as a quantity depending only on Cg , 77 and therefore ultimately on 
n,p, v, L. In turn we determine another integer Ki such that 

(3.14) 2"-f^i+2^^i < 1 , ^^2-^i''"+"-^i+2 < 1^ 

hold. This fixes Ki as a quantity depending on ri,Hi,Cd, and thcrcfore ultimately 
on n,p, V, L. 

Step 5: The degenerate case I. Wc consider the following situation: there exists a 
radius t < R such that p.9p happens to hold whenever r = t/2* and l < i < Hi G N, 
and finally we assume that both 

(3.15) \{Dv)Bt\ <'2VnM{2-"H) and s<2^M{2-"H) 

hold. We first notice that E{2-^H) < 2^M{2~"H), while iterating yields 

M{2~"H) < r]"^-Hl{t/2) . 
On the other hand we have that 

M{2-"H) < Tj"^-Hl{t/2) 



< CgT]^'-^ f {\Dv\ +s)dx 

JB, 



< Cgi^"''^ i \Dv-{Dv)BAdx + Cgif^-^[\{Dv)BA+s] 



< Cg-n^'-^Eit) + 4y^CgT]"'-^M{2-"H) , 

where we used p.lip . Therefore we have M{2^"H) < 2i^"^^^CgE{t), so that, again 
by (j3.13p . the final outcome is 

(3.16) E{2-"H) < \E{t) . 
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Step 6: The degenerate case II. Continuing the reasoning of the previous step, we 
assume that there exists a radius t < R such that p.9p happens to hold whenever 
r = t/2* and 1 < i < Hi + Ki N, assuming also, in alternative to (|3.15p . that at 
least one of the foUowing inequaHties holds: 

(3.17) \{Dv)bA > 2y/^M{2-"H) and s > 2y^M{2-"H) . 
In the case the first inequahty in (j3.17p holds we have that also 
\Dv - {Dv)bA > V^Mi2-"H) in B^-H.t ; 
thercfore, using also (|3.14p . we conclude as in (|3.16p : 

^(2-(^i+^i)i) < 2V^Mi2-^"'+'^'h) 
< 2y/^r^^^M{2-~"H) 



< 2if^j \Dv-{Dv)BAdx 
2nffi+i^Ki [ _ {Dv)b, I dx 

JBt 



< 

(3.18) < \E{t). 

On the other hand, if s > 2y^M{2-"H), we are in the situation of Step 3 with the 
choice r = 2~^^t and g = 2~(^i+^i)t; thercfore applying p.l2p yields 

E{2-^"^+^^h) < Cd2-^'P<'E{2-"H) 

while it also holds 

E{2-"H) < 2 / \Dv~ {Dv)b, I dx < 2""'+^E{t) , 

so that, merging the last two estimates and using (j3.14p wc again conclude as in 
(|XT5)) with 

Step 7: Summarizing Step 5 and Step 6. Here we first note that the result of 
Step 3 holds if we replace Hi by any larger integer since if (|3.13p holds for Hi , then 
so it does for any larger number. Thercfore, looking at Step 6; if (|3.9p holds for 
r = i/2* and 1 < « < i?i + Ki € N, and assuming also both 

\{Dv)bA < 2V^M{2-"'-^H) and s < 2y^M{2-"'-^H) , 

we can apply Step 5 with Hi replaced by Hi + A'i. Thercfore, summarizing the 
results of Steps 5 and 6 wc have that: If p.9p holds for r = i/2' and 1 < i < iJ := 
Hi + Ki then 

(3.19) E{Tt)<\E{t), t:^1/2". 

The cruciai fact is that H , and thercfore also t, only dcpcnds on i/, L. 

Step 8: Gonclusion. Lct use define (3i := 1/H so that (3i = f3i{n,p, ly, L) G 
(0,1] and r'^i = 1/2, where r € (0,1/2] has becn introduccd in (|3.19p : wc finally 
determino the exponcnt f3 appearing in (|3.6p by Ictting /3 :— niin{/3o, /?i} whcre Pq 
is the exponcnt from Lemma [321 Morcovcr, with p < iì as in the statement of the 
Theorem, wc fix e N such that t^^^R < g < t^R. Lct us now define 

S := {i e N : (EU) holds for r = i?/2' , i > 1} . 

We argue on an alternative; the first case is when § = N \ {0}. Then either Step 
5 or Step 6, and so Step 7, apply for every choice of the radius t = t^R and i > 0, 
thus obtaining that 

(3.20) E{t'R) < {l/2yE{R) = t^^^E{R) < t^^ E{R) 
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holds for every i G N. Taking into account the definition of /3i we have 
E{g) < 2t-''E{t^R) < ct^Pe{R) <c(^Y E{R) , 



so that (|3.6p foUows; note here that c = c(t) = c{n,p, v, L) by (|3.19p . In the other 
case we have that S 7^ (N\{0}), and therefore if we set m := min ((N\{0})\S) then 
((3^ imphes that |ì:»d| > M(i?/2""i) in the bah 5^/2™. At this point we observe 
that in the ball -8/2/2™ the matrix à{x) defined in (jS.lOp satisfles the foUowing 
uniform ellipticity bounds: 



c-i[M(iì/2'"-i) + s]f-2|A|2 < (a(a;)A,A), |à(2;)| < c[M(i?/2"-i) + s]p- 



2 



for c = c{n,p, v, L) > l and every choice of A £ R", and therefore we may apply 
Lemma 13.21 as already done in Step 3, thereby getting 



/3o 



(3.21) E{r)<ci^j^ j i?(i?/2") for every r < i?/2™ . 

Let iiow 7 € N be the unique non- negative integer such that < m < (7 + 
We apply Step 7 exactly 7 times thereby getting - also when 7 = 0- that 

(3.22) E{t'R) < {l/2yE{R) for every i < 7 . 
Moreover, using the definition of E{-) and t it immediately foUows that 

-B(i?/2™) < 2 -/ \Dv- ÌDv)b^-,J dx < 2'^"+^E{t^R) , 

SO that combining the last two inequalities gives 

E{R/2"') < T'"+^2-'< E{R) = ct"''^'E{R) < ct^^E{R) 

with c = c{n,p, V, L). On the other hand using ([^^ with r = r'f+^R = 2"(t+')^ì? 
and / G (N \ {0}), we gain 

< c2"I^°t'^'^°E{R/2"') 

< cT^^+'^l'EiR). 

The last estimate and p.22p in turn imply, for a suitable c = c(n,p,v, L), that 
the inequality E{t^R) < ct^^E{R) holds for every i G N. At this point we again 
conclude as after (j3.20p . □ 

3.3. Comparison estimates. This section is devotcd to the proof of a few suitable 
comparison estimates in which the density of the non-linear Wolff potential 

explicitly Comes into the play. We now fix, for the rest of Section 13.31 a ball 
B2R = B{xo, 2R) C fi with radius 2i?; we start defining w € u + Wq"^{B2r) as the 
unique solution to the homogeneous Dirichlet problem 



(3.23) 



div a{x, Dw) =0 in B2R 

w ^ u on dB2R . 



Lemma 3.3. Under the assumption (|1.2|) 9. let u G VF^'P(ri) he as in Theorem ì 1 . 1[ 
and w G u + Wq'^(ì?2/ì) as in (j3.23p . Then the following inequalities hold for a 
Constant c = c{n,p, v): 



(3.24) -i \Du~ Dw\dx <c 



WB2r) 
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(3.25) 



\u — w\ dx < c 



MB: 



2R) 



Proof. The proof revisit and modifies various comparison and truncation methods 
for measure data problems scattered in the hterature; a chief reference here is the 
work of Boccardo & Gallòuet [TOj. The proof is divided in three steps. We may 
without loss of generaUty assume that 2R = 1, that the bah in question is centered 
at the origin B2F(. = Bi, and finaUy that |/i|(i?i) = 1; the scaling technique necessary 
to reduce to such a situation is reported in Step 3 below. 

Step 1: The case p > n. We test the weak formulations of p.l|) and (|3.23P i 



(3.26) 



(a(x, Du) — a(x, Dw), Dip) dx 



ipdn , 



with if = u — w. Using Morrey-Sobolev's embedding theorem we estimate the 
resulting right hand side as foUows: 



{u — w) dfi 



< snp \u - w\[\n\iBi)] 

Bi 

I I 

< sup u — li; < 



\Du - Dw\P dx 



The last inequaUty used together with (|3.26|) with (p = u — w and (|3.2p yields 



\Du - Dw\Pdx < c 



Bi 



\Du~Dw\Pdx 



Bi 



so that we first get \\Du — Diì'IIìp(Bi) — ^^"^ eventually \\Du — Dw\\]^i(^g-^) < c 
that is ((XM)) - when = Bi and \fi\{Bi) = 1. 

Step 2: The case 2 < p < n. For fc > denoting an integer, we define the 
foUowing truncation operators 



(3.27) 



Tk{s) := max{-fc, min{s, k}} and <ì>fc(s) := Ti{s - Tk{s)) , 



defined for s G R. We test the weak formulation p.26p hy ip = Tk{u — w) . By 
setting Dk := {x ^ Bi : \u{x) — w{x)\ < fc}, we obtain, with c = c{n,p, v), 



(3.28) 



\Du - Dw\P dx < ck. 



D,, 



Similarly, testing p.26p with = $fc(w — u>) yields 



(3.29) 



/ \Du - Dw\P dx < c , 



where this time Ck {.t G -Bi : fc < \u{x) — w{x)\ < fc + 1} and c = c{n,p, v). By 
Hòlder's inequaUty, (|3.29p and the very definition of C^, we find 



\Du~Dw\dx < \Ck\~ 



\Du- Dw\Pdx 



Ck 



(3.30) 



p-1 n(p-l) 

< c\Ck\~ < ck 



-1 dx 



Ck 



We have of course used the elementary estimate 

1 



fc— 



\u — w] dx . 



Ck 
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Using (|3.30p with p.28p . Hòlder's inequality for sequences, and finally Sobolev's 
embedding theorem, with fco being a fìxed positive integer we bave 

f \Du~Dw\dx = f {■■■)dx+^ f {■■■)dx 

JBi JDko f.^^^ JCk 

fc-?(fcr) \u~w\'^ dxj 

y] / \u~w\^dx 



k—ko 



< ck^ + c 



(3.31) 



< cka+cH{ka) 



\Du — Dw\ dx 



In the last hnes we bave obviously set 

H{ko) := 

wbile c depends on n,p, v; note that 

P > 2 = 



E 

.k — kn 



n-jp - 1) 
n - 1 



> 1 



when fco — *■ oo- Now, ìi p < n 



so tbat H{ko) is always finite and satisfies i/(fc'o) 
then 

(P - 1)" ^ ^ . 
(n - l)p 

therefore we take fco = 1 in (|3.3ip and applying Young's inequality in (|3.3ip wc find 
(3.32) \\Du- DwWliìb,) <c, 

that is p.24p when 2i? = 1 and |/.i|(i3i) = 1. When p = nwe choose fco = ko{n,p, v) 
large enougb in order te bave cH{ko) = 1/2 in (|3.3ip and reabsorb the last integrai 
on the right band side, so that p.32p foUows again. 

Step 3: Scaling procedures. We first reduce to the case B2R = i?i by a standard 
scaling argument. i.e. letting 

u{xo + 2Ry) w{xo + 2Ry) 



ù{y) 



2R 



2R 



and 



à(y, z) := a{xo + 2Ry, z) , jì[y) 2Rfi{xo + 2Ry) 

for 2/ e i?i so that — div à(y, Du) ~ fi and div à{y, Dw) ~ bold. At tbis point 
one writes estimate p.24p for ù and w and then scales back. To reduce to the case 
1^1 (Si) = 1 we adopt another scaling: tbis timo we define A :— [|/ì|(ì?i)]^/(p~^\ 
and we may assume ^ > otherwise u = w and (j3.24p foUows trivially by the 
strict monotonicity of the operator a(-). We define the new solutions u := A^'^u, 
w := A~^w, the new datum /2 := A^~Pp,, and the new vector field d{x, z) := 
A^~Pa{x, Az). Therefore we bave div d{x, Du) = /2, div d{x, Dw) = 0, in the weak 
sense. We make sure that we can apply the result in Step 2. Trivially \fl\{Bi) = 1 
and morcover it is easy to see that the vector field d{x, z) satisfies (ll.2p with s 
replaced hy s/A>Q. Therefore p.24p bolds in the form 



\Du — Dw \ dx < C2 
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with C2 = C2{n,p, v). At this point we find back (j3.24p just using the definitions of 
u, w. As for (|3.25p . this is a consequence of (j3.24p via the use of Sobolev inequahty: 



+ \u — wl dx < cR -t \Du — Dw\ dx < c 
j B2R ^ B2R 



2R) 



R"-P 

The proof is now complete □ 



Remark 3.1. In the proof of the gradient estimate p.6p we shall only need (j3.24p . 
while (|3.25|) will bc nceded only later, in the diffcrcnt contcxt of Remark |4 . 1 1 bclow . 

After introducing w in p.23p wc similarly define v ^ w + Wq'^{Bjì), on the 
concentric smaller ball Bjì = B{xo, R), as the unique solution to the homogeneous 
Dirichlct problem with frozcn coefficients 

div a(xo, Dv) = in Bn 

V = w on dB /f . 



(3.33) 

Next, we state a standard comparison result; we report its proof for completeness. 



Lemma 3.4. Under the assumptions of Theorem \l.ll with w as in (|3.23p and v as 
in (j3.33p . there exists a Constant c = c{n,p,v,L) such that the following inequality 
holds: 

(3.34) / \Dw - Dv\P dx < c[LiLo{R)Y i {\Dw\+sYdx 

J Br ^ Br 



Proof. Wc note that by the growth condition on a in (|1.2[) i and by p.2p . we also 
have the following coercivity condition: 

(3.35) c-^{\z\^ + s^)'^\z\^ - csP < {a{x,z),z), c = c{n,p,v, L) > l , 

which easily foUows by Young's inequality. Using (|3.35p . that is testing (|3.33P i with 
w — V, and thcn (|1.2p 9 and Young's inequality, gives the energy bound 



(3.36) 



f \Dv\Pdx <c ( {\Dw\ + sYdx, 

J Br Br 



which holds for a Constant c depending on n,p, i^, L. In turn, using (j3.ip . the fact 
that both V and w are solutions, (|1.2p ^ and again Young's inequality. we have 

( {s^ + \Dv\^ + \Dw\^)^\Dw - Dv\^ dx 
Jbr 

< c I (a(.To, Dw) — a{xo, Dv), Dw — Dv) dx 

i 

{a{xo, Dw) — a{x, Dw), Dw — Dv) dx 

t 

< cLiLj{R) I {\Dv\^ + \Dw\^ + s^)'^\Dw - Dv\ dx 
Jbr 

<\ f {\Dv\^ + \Dw\^ +s'^)'^\Dw- Dv\^dx 
2 Jbr 

+ c[LiUj{R)f [ {\Dv\ + \Dw\ + s)P dx . 



Therefore we gain 

/ {s^ + \Dv\^ + \Dw\^)'^\Dw~Dv\^dx <c[LiUj{R)]^ f {\Dv\ + \Dw\+ s)p dx, 

JB^ J Br 

and (|3.34p foUows by last inequality rccalling that p > 2. □ 
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Lemma 3.5. Let u he as in Theorem and let w ^ u + Wq'^{B2b) md v £ 
w + Wq''^{Br) be defined in (|3.23p and (|3.33p . respectively. Then for a Constant c 
depending only on n^p,v^L, it holds that 



(3.37) 



i \Du-Dv\dx < c\l + [LiUj{R)]p] 
Jbr ^ 

+c[LiUj{R)]ì 



\Du\ + s)dx . 



Proof. We start proving the foUowing inequality: 



(3.38) 



Dw- Dv\dx < c[LiUj{R)Yi 



{\Dw\ +s)dx, 



B2R 



with c = c{n,p,v, L). Keeping (j3.35p in mind we may apply Gehring's lemma 
in the version presented in [MI Chapter 6], finding there exists a Constant xo = 
Xoii^yP-, ì^, L) > 1 such that the reverse Hòlder type inequahty 

(/ i\Dw\ + s)''''P dx] <ci {\Dw\+sYdx 

holds whenever Bg C B2R, for a Constant c depending only on L. In turn, 

applying Lemma [XT] with g = {\Dw\ + s)^, leads to establish that also 



(3.39) 



{\Dw\+sYdx\ <c 



B^R 



D'w\ + s) dx 



holds. Using now p.34p and previous inequality we estimate as foUows: 



\Dw — Dv\ dx < c 



\Dw- DvFdx 



< c[LiLu{R)]p i^j- {\Dw\+s)Pdxj <c[LiUJÌR)]p 
and therefore (|3.38p foUows. Using p.38p and (|3.24p now yields 

MB2B 



\Dw\ + s)dx 



B2R 



(3.40) 



\Du — Dv \ dx < c 



R" 



+ c[LiUj{R)p4 {\Dw\+s)dx. 



In order to estimate the last integrai in (|3.40p we simply use p.24p as follows: 



\Dw\dx < 



< 



j- \Du\dx + l 
J B2R ^ B2R 

\Du\ dx + c{n,p, u) 



Du — Dw\ dx 

t^\{B2R) 



B2R 



R" 



Using the last inequality in combination with (j3.40p we conclude with p.37p . □ 

3.4. Proof of the main estimates. We start with a technical lemma. 

Lemma 3.6. Let u he as in Theorem ll.ll then there exists constants (3 G (0, 1] and 
c, Ci > 1, ali depending only on n^PjU^L, and a positive radius Ri = Ri{Li,uj{-)), 
such that the following estimate holds whenever Bg C Bj^ C B2B C SI are concentric 
halls with R < Ri: 

/ \Du-{Du)B^\dx<ci(%Y -f \Du-{Du)B2R\dx 

Jb, JB2R 



(3.41) 



R 

-c I — 



Du\ + s) dx . 
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Moreover, in the case the vector field a(-) is independent of the variable x, the 
previous inequality holds without any restriction on R. 

Proof. Starting by B2R we define the comparisoii functions v and w as in (j3.33p 
and (|3.23p . respectively, and then we compare Du and Dv by mean of (|3.37p . using 
as basic reference estimate for v, eventually transferred te u: 



< c 



/ \Du-{Du)B^\dx<2-f \Du-{Dv)B^\dx 

< 2 4 \Dv~ {Dv)b^ \dx + 2-l \Du- Dv\dx 
Jb^ " JBg 

<cf-|) j \Dv~iDv)Bn\dx + ci^-j j \Du-Dv\dx 
S;Y -f \Du-(Du)BJdx + c(-\ i \Du-Dv\dx 

JBn Ve/ JBn 

(3.42) <2''+^c(^Y -f \Du- {Du)B,Jdx + c(-] / \Du~Dv\dx. 

Notice tliat from tlic second-last to the fast line we estimatcd as follows: 
i \Du- {Du)Ba\dx < i \Du- {Du)B2n\dx 

(3.43) +\{Du)Bn-iDu)B,n\<2''^'-f \Du^ iDu)B,n\dx. 

J B2R 

In order to gct (j3.4ip it is now sufficicnt to estimate the fast integrai in (j3.42p by 
mean of (|3.37p and then to take the radius Ri such that LiU!{Ri) < 1. Necdlcss to 
say, in the case a(-) does not depend on the variable x we can take lu{-) = and 
therefore no restriction on R is needed. □ 

Proof of Theorem \l.l\ The proof is dividcd in three steps. In what follows ali the 
radii considered will be smaller than a certain radius R: 

(3.44) R<R 

whose final values will be determined towards the end of the proof, i.e. we shall 
dccrease the valucs of R scveral times according to our nceds, but always in 
such a way that the rcsulting detcrmination of R will be stili depcnding only on 
n^Pyiy, L, Li,uj{-). This will finally givo the radius R in the statement of Theorcm 
11.11 From the proof it will be clcar that such a restriction is ncccssary only in the 
case the vector ficld a(-) dcpends on the variable x, while ali the incqualitics holds 
with no restriction on R when a = a{z); the final outcome is that estimate (|1.6p 
holds for cvery ball B{xo, 2R) C fi. This observation also clarifies the last assertion 
in the statement of Theorem 11.11 Initially, we shall take R < Ri whcrc the radius 
Ri = i?i(Li,a;(-)) has been determined in Lemma 13.61 abovc. In the following we 
shall consider sevcral constants, in general depcnding at least on the parameters 
n,p, L; rclevant depcndence on additional parameters will be emphasizcd. 

Step 1: Basic dyadic sequence. Referring to estimate (|3.4ip . we select an integer 
H = H{n, p, VjL) > 1 such that 

1\^ 1 



(3.45) ci^-) <-. 

We notice that the depcndence of H upon n^p,v^L comes from the similar de- 
pcndence of (3 and ci presented in Lemma 13.61 Applying (j3.4ip on arbitrary balls 
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Bg — Bji/2H ^ ^_R/2 C Br and using the fact that uj{-) is iion-decreasing we gain 
/ \Du^{Du)B,,,„\dx<-i \Du^{Du)BAdx 



(3.46) 



- C2 [LiUJ 



(I 



Du\ + s) dx , 



where c, C2 depends only on L, H and thcrcforc ultimatcly on i/, L. By 

the elementary estimation 

i {\Du\ + s) dx < j- \Du-{Du)Bjdx+\{Du)Bj + s 
estimate p.46p turns to 

1 



R/2H 



Du-{Du)BR,,Jdx< [-+C2[LMR)]^ f \Du - {Du)BR\dx 



(3.47) 



4 



c[LMR)]-p {\{Du)br\ + ^ 



We start reducing the value of R: we take R, depending only on n,p,v, L, Li and 
w(-), sniall enough in order to get 



(3.48) 



Notice that in order to estabUsh the claimed dependcnce of R upon the various 
parameters we have used that C2 depends in turn on v, L. By (j3.44[l . as uj{-) is 
non-decrcasing. mcrging (|3.47p and p.48p yields 



(3.49) 



/ \Du^{Du)BR,,„\dx<]-j \Du-{Du)BR\dx 

JBri-iu JBr 



+c 



+ c[LMR)] - {\{Du)br\ + s) 



We now fix a bah B{xo,2R) C fi as in the statement of Theorcm 11.11 for i G 
{0, 1,2, . . .}, let US define 

(3.50) B, := B{xa,R/{2Hy) := B{xo,Ri) and h \{Du)b^. 
For every integer m G N we write 

fc,„+i = ^(fc,;+l " fci) + fco < ^+ iDu - (Du)bJ + fco 

i = i = Q''B^ + l 

(3.51) < V(2iJ)"-/ (i:'u)sjda; + /c( 



• 



i=0 



Therefore, defining 

(3.52) Ai -.^l \Du- {Du)B^dx 

J B, 

inequaUty p.5ip becomes 



(3.53) 



< (2i/)" ^ A, + fco . 
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To estimate the right hand side of the previous inequahty we observe that (|3.49p 
used with R = yields, whenever i > 1 



(3.54) A,<^A,^i+c 



n — l 
i-1 



+ c [LiUj{Ri^i)] p (fc^_i + s) . 



with c = c{n,p, ly, L), where we have taken into account (j3.50|l and that H dcpcnds 
on n,p, v, L. We now consider (j3.54p for «€{!,..., m} and sum up over i, thcreby 
gaining 



m— 1 



m — 1 



i=l i=0 1=0 

and therefore 

m m — 1 



[LMR,)]Ì {h + s) , 



1=0 



i=0 



i=0 



Using the last inequahty in (j3.53p yields, for every integer m > 1 



m — 1 



(3.55) fc™+i <cÌAQ + ko+J2 



i=0 



in — 1 



- c^[LiUj{R,)]p {k, + s) , 



i=0 



as usuai for a Constant c = c{n,p, ly, L). 

Step 2: Wolff Potentials and Dini-continuity. We notice that 

m — 1 ri 



E 



RI' 



< 



< 



< 



oc 

E 


mB:) 


[ Rr' 


i=0 




■ p2R - 


iog 2 Jj, [ 


2^ 


■ p2R - 


log2 7fl [ 



{2H)-. 



\li\iB{xo,g)) 



\li\iB{xo,g)) 

gn-l 



E 



log2i? — 



"-^ dg v-^ 
— + > 

^ 1=0 

1 

"-^ dg 

g 

'\p,\{B{xo,g)y 



g 



,11-1 



'i+l 



-1 dg 

g 



< 



2^ 



Ì2H)-- 



Wt {xo,2R). 



log 2 log2ff 

Rccalling the dependence of H determined in p.45p we conclude that there exists 
a Constant c depending only on 7i,p, v, L such that 

m— 1 

(3.56) 



4=0 



WBi 



<cW'i (xo,2R) 



holds whenever m € N. As for the last tcrm in (|3.55p . using the fact that uj{-) is 
non-decreasing, we estimate 

m — 1 oo 



i=0 



< 



< 



1=0 

1 

fog2 
1 

fog2 



2R 



Hg)]i^ + J2HRi+i) 
g ^ 



2R 



>{q)V' 



i=0 

dg 1 



g \og2H 



E 

i=0 



dg 



H9)] 
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Recalling the definition of d{-) in (|1.3p . and the fact that > 1 we have 

(3.57) •£k«,,|Ì<MM. 

For later use, we now further restrict the vaine of i? in order to have 

(3.58) Lld{2k)<\. 

In light of incquahtics (|3.56|) and (j3.57p . (|3.55|) yields, for every m > 1 

m— 1 

(3.59) kra+x < C3M + C4 ^ [LMR,)]h, , 

i=0 

where C3, C4 depend only on n,p, v, L and we have set 

(3.60) M:=/ {\Du\ + s)dx + W't {xo,2R) . 

Jb{xo,R} "'^ 

Notice that, using p.58p and Hòlder's inequahty, to get p.59p from (|3.55p we have 
also estimated 

(3.61) Aa + ko + d{2R)s<c4- {\Du\ + s)dx. 

Jb(xoM) 

Finally, as in (|3.5ip . by estimating fci = {ki — /cq) + < cAq + fco, and using l|3.6ip . 
we complement p.59p with 

(3.62) fco + fci < C3M 

that we obtain by enlarging C3 a bit. 

Step 3: Induction and conclusion. We restrict the value of R for the last time in 
order to have 

(3.63) d{2R) < ^ 



where C4 = Ci{n,p^ v, L) is the Constant appcaring in p.59p : this. togcther with 
the choices p.48p and (j3.63p . finally deterniines i? as a positive quantity depending 
only on 7i,p, i^, L, Li, a;(-), as required in the statement of the Theorem. We now 
prove that the foUowing inequality holds whenever m > l: 

(3.64) fc,„+i<2c3M, 

where C3 and M have been introduced in (j3.59p and (j3.60p . respectivcly. We prove 
p.64p by induction; the cases m = —1,0 are settled in (|3.62p : now we assume the 
validity of (|3.64p whenever m < m, and prove it for m + 1. We have, using (|3.59p . 
(12111), (ESZl), and 

m 

^rh-\-2 — Y,[Liu;iR,)]Ìh 
i=0 



< 



C3M + 2C4C3M ^[Liw(iìi)]' 



< c3M+^^^^^di2R)M 

log 2 

< C3M + 8c4C3Ll^^d{2R)M 
(3.65) < C3M + C3M = 2C3M , 
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SO that (|3.64p holds whenever m £ N. We now recali that, being Du assumed to 
be a eoiitinuous veetor field, for every xo G SI it holds that 

(3.66) \Du{xa)\^ lim fc,„+i < 2C3A/ . 

m — >oo 

Merging the last inequality with (|3.60p finally yields (jl.6p . For the easc with no 
x-dependence in the veetor field a(-) we just notice that in this case the last term 
in p.46p does not appear, and the same happens in (|3.49p . therefore no smallness 
assumption of the type p.48p or (|3.63p are required, and ali the estimates in the 
iteration procedure work with no restriction of the type (|3.44p . The final outcome 
is that estimate (|1.6p holds for every ball B{xq,2R) C fi. We finally observe that 
the only point in this proof where we used that w e (fi) was at the very end, 
in order to assert (|3.66p at every point xq; otherwise we could have just assumed 
u € Wi^^{iì) to get p.66p almost cverywhere. □ 



4. Parabolic estimates 

This section is devoted to the proof of Theorems 11.21 and 11.41 and therefore 
we shall in general argue under assumptions (jl.lSp and p.28p . respectively. We 
observe that assumptions (|1.18P i- (jl.l8p 9 obviously imply (|1.28p . and this fact will 
be implicitly used several timcs in the foUowing. 

4.1. Basic estimates from Nash-Moser's theory. Here we shall emphasize a 
few properties of solutions ù e C'^(ii, Ì2; i^(B(x, 7))) n i^(ii, <2; W^^'^(5(x, 7))) for 
7 > 0, to homogeneous , non-linear, parabolic equations of the type 

(4.1) ùt - div ò(a;,t,L»M) 0, 

therefore considcred in the basic cylinder Q = B{x,j) x (ti,t2)- The veetor field 
b: Q X M" M" is supposed to be only Carathcodory regular - in particular 
the dependence on the coefficients (x, t) is mcrcly mcasurablc. The ncxt rcsult 
essentially cncodcs Nash-Moser's regularity for parabolic equations. 

Proposition 4.1. Let ù e C"{ti,t2; L'^{B{x,"f))) n L'aiti, t2;W'^''^{B{x,-/))) be a 
weak solution to (j4.ip . under the assumptions (jl.28p . Then u S Ciof (Q) for some 
(3 G (0, 1] depending only on n, v, L, and moreover, there exists a Constant c, again 
depending only on n, v, L, such that the following inequality: 



(4.2) / \ù-{ù)Q^\dxdt<c(j-)^i 

jQe JQr 



\ù — {u)qj^ \dxdt + csR 



holds whenever Qg C C Q are parabolic cylinders with the same vertex . 

Proof. We shall rely on the precise estimates given by Lieberman in |411 Chapter 
6]; this is in turn an extension of the classical Nash-Moser's theory for linear par- 
abolic equations [51] [52] . We first observe that by a simple application of Young's 
inequality from (|1.28p it also follows 

(4.3) (i^/2)NP < {b{x,t,z),z) +€3"^, 

where c = 2L'^ /v; compare with p.35p . Therefore applying [IT] we have 



(4.4) ose ù < c — — ) ose ù -1- csR for every g < R/2 , 

Qg \R/2 J Qr/2 

for constants c > 1 and (3 G (0, 1], depending only on n, L. Here we are adopting 
the standard notation 

ose ù = sup |ù(a;, t) — ù{x, t)\ {x, i), (.t, t) e Qg . 

Qe Q„ 
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Estimate (|4.4p maybe neeessitates a few elarifieations; it can be inferred from [JTJ 
Theorem 6.28] once the following things are taken into account. First, one has to 
adapt the estimates derived in [411 Chapter 6] , considering the growth assumption 
(|4.3|) and to compare this with [4TJ (6.20)] and j4ll (6.28)]. Second, the estimates 
in [m Chapter 6] are often stated for hnear operators, but the hnearity is actually 
irrelevant, as growth conditions are the only relevant thing, therefore aU the esti- 
mates can be derived for general quasihnear operators; see the comments at the 
beginning of [U Section 6.5]. Again from [?TJ Theorem 6.17] wc gain 

(4.5) sup 1m1<c4 \ù\dxdt + csR, 

Qr/2 JQr 

again for a Constant c depending only on n, ly, L, and whenever Qn C Q. With 
m £ M being a fixed number, we observe that ù — to stili solves (|4.ip . therefore 
estimate (|4.5p gives 



(4.6) sup ]ù — m]<c4 

Qr/2 C 



\ii — to] dx dt + csR . 

Qr/2 JQr 
On the other hand wc have 

(4.7) ose (ù — to) < 2 sup ]ù — to] < c -f ]{t — to] dx dt + csR , 

'^«/a Qr/2 Jqr 

and. in turn applying estimate (|4.4p to ù — to wc find, by mean of (|4.7p 

ose ({t — to) < c f-^^ -[ \ù — m\dx dt + csR . 
Q. Jq„ 



ose u 

Qe 



Noticing that 



we finally obtain 



/ 

JQ. 



\v, — (u)q^ \ dxdt < ose ù 



~f \ù~{u)q \dxdt<c(—^ -f \ù — m\ dx dt + csR . 

Jq, ' Jqr 

At this stage (j4.2p foUows choosing to ~ {u)qj^ . □ 



4.2. Comparison estimates. In the rest of Scction [4?2l we keep fixed a symmetric 
parabolic cylinder Q{xo,tQ,2R) = Q2R C ÌÌt- Wc start considering the uniquc 
solution 

(4.8) w e C"(to - 4i^^ to; L^{B{xo,2R))) n L^{to - 4i?2, to; W^'^{B{xo,2R))) 
to the following Cauchy-Dirichlet problem: 

J wt ~ div a{x,t, Dw) ~ in (52_r 

I W = W on dpaiQ2R ■ 



(4.9) 



Lemma 4.1. Let u e C°(-r, 0; L^{VL))C]L^{-T, 0; W^^^{n)) he a solution to (fTTTp 
with /i G L'^(Q,t), under the assumptions (|1.28p . and /ef w be defined in (|4.8 
Then there exists a Constant c ~ c{n, v) such that the following estimates hold: 



{ 

JQ2R 



(4.10) f \Du- Dw\dxdt< 



c\ì^\{Q2r) 

JlN-l 

(4.11) i \u~w\dxdt < 

JQ2R 



c \t^\iQ2R) 
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Proof. Step 1: Universal estimate. Here we assume {xo,to) = (0,0) and 2R = 1 - 
that is Q2fl = Qi - and that |/i|((5i) = 1, and prove that the foUowing universal 
inequality holds: 

c{n) 



(4.12) 



\Du- Dwìdxdt < 



Our starting point here wiU be the parabohc estimates developed in [9]. As de- 
scribed in Section 2, we use the Steklov-averages formulation of both (|1.17p and 
(|4.9P i ■ i.e., for every t e (-1, 0) and t + h<0 there holds 



(4.13) 

and 
(4.14) 



(j3tUhi-,t)ip + {[a{-,t,Du)]h,Dip)^ dx = / ipfih{x,t)dx 



{dtWh{-,t)ip + {[a{-,t,Dw)]h,Dip)^ dx 







for every ip G C^{B) with compact support, and by density whenever Lp G Wq''^{B). 
The initial datum of Wh is here taken in the sense of L'^{B), which means that 
Wh{-,-l) u{-,-l) in L^{B) when /i i 0, and since u e C°([-l, 0); ^^(B)), this 
implies in particular that 



(4.15) 



l|m||(u/i - u;/0(-,-1)I1l2(b) . 



With the notation fixed in p.27p . by ^'fe : R ^ M we denote the following primitive 
'^k{s) ■= /q $fc(K) dn for r e M, that for later use we compute explicitly: 



' i + (r-fc-l) 
i(r-fc)2 

*fc(T) = l 

i(r + fc)2 
i - (T + fc + 1) r < -fc- 1 



T>k + l 
k < T < k + 1 
-k <T <k 
-k - 1 < T < -k 



and note that "i/k > 0. By testing the difference Steklov-formulations (|4.13p . (|4.14p 
J (dt{uh - Wh){x, t)ip + {[a{-,t, Du)]h - [a{-,t, Dw)]h,Dp)^ dx 

(4.16) = / ippih{x,t)dx , 

Jb 

with the choice f{x,t) :— C{t)^k{uh ~ 'Wh){x,t), x E B, where C(') is a smooth 
function, and then integrating the resulting equality on (—1,0) with respect to t 
we obtain 

dt{uh-Wh)^kiUh-Wh)C dx dt 

{[a[-,Du)]h - [a[-,Dw)]h,D^k{uh-Wh))C dx dt 

^k{uh-'Wh)CtJ'h dxdt. 

i) as foUows: 



(4.17) 



With re (-1,0) now choose C e 



1 K < T 

Ce{K) = CÌK)={ l-^{t-K) T<K<T + e 

K> T + e . 
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Then, the first integrai in (j4.17p can be rewritten in the forni 
dtiUh-Wh)^k{uh-Wh)C dx dt 

dt[i>k{uh-Wh)C] dxdt - / k{uh-Wh)C,t dx dt 
' JQ 

(4.18) =-/ ^k{uh~Wh){x,-l)dx- k{uh-Wh)C,t dx dt . 



Taking into account the special choicc of C wc sec that the second integrai appearing 
on the right-hand side of (|4.18p converges to /g ^k{uh—Wh){x, r) dx when e | for 
almost every r € (—1,0), whilc the first integrai vanishes as ft. | in light of (j4.15p . 
As we have seen above the initial datum is taken in the sense of L^{B). Combining 
with ((ÌTT)) and letting first e [Q and then /i i in wc obtain 

^k(u~w)(x, t) dx 

{a{x,t, Du)—a{x,t, Dw), D^kiu—w)) dx dt = / / ^k{u~w)iJLdx dt 
iJB J-iJb 

for almost every every choice of r e (—1,0). This Icads to 



sup / kiu--w){x,T) dx 

-1<t<oJb 



{a{x, t, Du)—a{x, t, Dw), D^kiu~w)) dx dt 

Q 

(4.19) < / \^k{u-w)\\^l\dxdt. 



We first exploit the properties of 'i!k in order to estimate the first tcrm of the 
left-hand side from bclow. We start with the obvious estimate 



(4.20) sup / *o(M-w)(-,r)da; < l/^KQ) < 1, 

-1<t<0Jb 

which foUows by applying (|4.19p for k = 0; note also that |m|(Q) = 1 by assumption. 
With T being fixed, decomposing 

B^{xe B : \{u-w){x,t)\ < 1} U {x e B : \{u~w){x,t)\ > 1} =: B'^UB'^ 

we have for almost every t e (—1,0) that 

'^o{u—w){x, t) dx 

ì 

(4.21) = / ì;\{u-w){x,T)\''dx+ [ \{u~w){x,T)\dx- ^\B';\ 

J B'^ " JB!; 

holds. On the other hand we trivially have /g, \{u—w){x,T)\dx < \B\ , so that 
merging the last estimate with (I4.20|) - (|4.21[) we have 

\{u-w)ix,T)\dx < + |/i|(g) < c(n) 

which holds again for almost every t e (—1,0), and thereforc wc conclude with the 
L°°-Li-bound 

(4.22) ||w - ■u;|lic»(_i.o:Li(s)) sup / \{u - w){x,T)\dx < c{n) . 

-1<t<0Jb 
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We now start exploiting the second integrai appearing on the left-hand side of 
(gm. For this we define 

Ak := {{x, t) eQ : k< \u{x, t) ~ w{x, t)\ < k + 1} . 

Then, using the fact that D^k{u — w) = Du — Dw on Ak and D^k{u — w) = 
otherwise together with (|1.28p 9. $fc(-) < 1 and \^J.\iQ) = 1 in (|4.19p we obtain 

\Du - Dw]"^ dx dt < {a{x,t,Du)-a{x,t,Dw),D{u-w))dxdt <1 
so that the estimate 

(4.23) / \Du - Dw\^ dx dt < - 

holds for every fc g {0, 1,2,.. .}. Now, for A > 1 we estimate 
12 ^ r \Bu-Dw? 



f \Du-Dw\'^ , , ^ f L 
/ T \ '--^dxdt=y T- , 



dx dt 



k=0 ' 

oo 



(4.24) < g /^^ \Du - Dn^f dxdt < . 

Here we have used in the last Une the uniform estimate (|4.23p . This allows us to 
argue as foUows: For ahiiost every r G (—1,0) we have by the Cauchy-Schwartz 
inequahty that 

< f / '''ir-'^^rl. ^-^ìV / icin) + Ku^.)ix,r)\)^dx ^ 



{1 + \Ìu-w){x,tW 
Integrating the preceding inequahty with respect to t on (—1, 0) and using Cauchy- 
Schwartz with respect to r together with (|4.24p then yields 

\Du~Dw\ dx dt 



^ ( / J^^y ""^1" dxdtVf f I {^ + \{u~w){x,t)\fdxdt 



(4.25) < 



{l + \{u-w){x,tW 
(1 + \u-w\)^ 

^\\{u-w){xMÌ>-(B)dt) 



Q 

2^-iA r / rO 



KA-1) 



At this stage we use a wch-known version of the Gaghardo-Nirenbcrg embedding 
thcorcm which rcads in our sctting as foUows: 

\\{u~w){x,t)\\LX^B) < c(?^)||£'(^t-w)(2;,^)||^l(B)||(^^-w)(x,^)||^7(B) , 
which holds for every choice of 9 and A satisfying 

1 9 

0<9 <1, l<A<oo, - = 1- -. 

A n 

We note that we have A > 1 at our disposai to ensure 9 G (0, 1). Recalling the 
L°°— L-'^-estimate (j4.22[) the second integrai on the right-hand side of the preceding 
inequahty is bounded by c(n)^~^, and therefore 

r \\{u~w)ix,t)\\l.^s^dt < c{ny+^^'-'^ f \\D{u-w){xML^B)dt. 
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We now perforili the choice of and A. In order to have 6\ = 1, i.e. 6' = j, the 
identity relating 9 and A yields A = so that 9 = This imphes that 

\\(u-vj){x,t)\\^\\ dt<c{n) \D{u-w)\dxdt. 

^1 L >. (S) Jq 



Inserting this last estimate in (j4.25p , and using Young's inequahty we arrive at the 
foUowing inequahty, which is turn imphes (|4.12p 



\Du- Dwìdxdt < 



< 



1 



c(n) 
c{n) 1 



Du — Dw\ dx dt 
Du — Dw\ dx dt . 



Step 2: General case and scaling. Similarly to the eUiptic case we first reduce 
to the case Q2R = QiXo,to;2R) = Q hy changing variablcs and passing to new 
Solutions and vector ficlds. More precisely, for {x,t) S Q 



(4.26) 



ù(i, t) 
w^x, t) 



u{xo + 2Rì,ta + mH) 
2R 

w{xo + 2Rx, to + AR'^i) 
2R 



à(i, t, z) := a{xo + 2Rx, io + 4i?^<, z) 
. t) 2Rn{xo + 2RÌ, to + AR^t). 

Then it follows that — div à{x, t, Du) = jl and wj. — div à(i, t, Dw) = in Q, and 
ù = ù; on dparQ- Furthermore the new vector field à(-) satisfies (|1.28p . To remove 
the additional assumption l/^KQ) = 1 we assume, as we can, that A := \ii\{Q) > - 
otherwise there is nothing to prove since the strict monotonicity of the vector field 



would imply u 



and then we re-scale u, w, a(-) and fj, as follows: u := A 



w := A~^w, a(x, t, z) := A~^a{x, t, Az), p, := A~^^. Then, it is easily verified that 
M is a solution of ut — diva(x, t,u) = fi on Q, that w solves wt — diva(a;, t,w) = 
on Q, that u = w on QparQ and l/^KQ) = 1. Moreover, a(-) fulfills the strict 
monotonicity condition (jl.28p which is the only onc uscd in a quantitative way 
in the dcrivation of the univcrsal comparison estimate (j4.12p . Thcrcforc estimate 
(j4.12p applied to u ~ w yields 



/ \Du - Dw\ dx dt < c{n) 
Jq 



iQ 

and re-scaling back this inequahty from u — Hi io u — w Icads us to the compar- 
ison estimate (j4.10p with R = 1, and by the previous scaling (|4.10p is proved 
in the general case. In order to get (|4.1ip it is sufScient to recali that {u — 
w){-,t) £ Wo' {B{xo,2R)) for almost every t G {to — 4i?^,to) and therefore ap- 
plying Poincaré's inequahty slicewise and integrating we gain 

/ \u — w\dxdt<cR / \Du — Dw\dxdt 

JQ{xo,to;2R) J Q{xo,to;2R) 



lQ{xo,to;2R) 

which together with (|4.10p implies (|4.1ip . This complctcs the proof. 



□ 



We now proceed with a further comparison estimate; after having defined the 
comparison solution w in (j4.8p we define the unique solution 



(4.27) V G C°(fo - R^ to; L\B{xo.R))) n L^to - R^to; W^'^{B{xo, R))) 
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of the foUowing Cauchy-Dirichlet problem: 

J Vt ~ div a{xo,t, Dv) ^ ìt^Qr 

\ V On ^parQiì ■ 



(4.28) 



We recali the reader's attention on the fact that in (|4.28|) we have frozen the coef- 
ficicnts only with respect to the space variable. 

Lemma 4.2. Lei u e C"(-r, 0; L'^{n))nL^{-T, 0; W^''^{n)) he a solution to p:T7|) 
with fi G L^{ÌÌt), under the assumptions (|1.28p . let w he defined in (|4.8p - (|4.9p . and 
finally let v he defined in (|4.27p - (|4.28p . Then there exists a Constant c = c{n,v,L) 
such that the following estimate holds: 



f \Du~ Dv\dxdt < c[l + LiUj{R)] 

Jqr 



lMl(Q2fl) 



iQr 

(4.29) +cLiLu{R)f {\Du\ + s) dx dt . 
Proof. Let us first get the estimate 

(4.30) / \Dv - Dw\'^ dx dt < cLl[u}{R)Y 4 {\Dw\ + sf dx dt . 
Jqr JQr 

The proof consists of a rathcr standard comparison argumcnt which we report here 
for the sake of completeness. The following computations are formai, as they would 
need the existence of the time derivatives for both v and w; on the other hand they 
can be made rigorous by using the Steklov-averages formulation exactly as done for 
(|4.17p - (|4.18p . Using the fact that both v and w are solutions yields 

{v — w)t — div {a{xQ, t, Dv) — «(xo, t, Dw)) 

(4.31) = div (a(xo, t, Dw) — a{x, t, Dw)) . 

By testing the weak form of the previous equality with v — w - here we need to pass 
to Steklov averages - we get, after standard parabolic manipulations and the fact 
that V and w agree on 9pai QiÌ! that 

sup / \v — w]'^ {x , t) dx 

to-R?<t<to Jb{xo,R) 

{a(xo, t, Dv) — a{xQ, t, Dw), Dv — Dw) dx dt 

{a{xo, t, Dw) — a{x, t, Dw), Dv — Dw) dx dt 



< 

JQr 

Discarding the first term in the previous inequality, using (|1.28p and p.lSP '^. yields 

i \Dv - Dw\^ dx dt < cLiUj{R) 4 {\Dw\ + s)\Dv ~ Dw\dx dt . 

JQr JQr 

At this stage (j4.30p foUows by the last estimate and a standard use of Young's 
inequality. The next step is to recali a higher integrability estimate which foUows 
form the application of Gehring's lemma in the parabolic setting. The statement 
can be found in several paperssee for instance |53| . once taking into account the 
growth conditions p.lSp considered here. There exist constants c > 1 and x > 2, 
depending only on n, v, L, such that the following inequality: 

Ij \Dw\^dxdt\ <c\J {\Dw\ + s)"^ dx dt\ 
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holds whenever Q2g C Q2R- We are tlierefore in position to apply Lemma 13.11 
which allows to establish 

(4.32) (/ \Dw\^ dx dt] < c 4 (\Dw\ + s) dx dt . 

\JQr ) JQlR 

Combining (|4.30p and (|4.32p via Hòlder's inequality we gain 

4 \Dv - Dw\dxdt < LiUj{R) 4 {\Dw\ + s) dx dt . 

JQr JQ2R 

Again combining the last estimate with (|4.10p . and actually using it twice, finally 
gives (11211). □ 

We conclude this section with yet another a priori estimate, in which, for the 
first and only time in this section we use the complete assumptions (|1.18p rather 
than the wcakcr oncs in p.28p . 

Proposition 4.2. Let v he defined in (|4.27p - (|4.28p .' there exist constants c > 1 and 
(3 € (0, 1], depending only on n, v, L, such that the following inequality: 



(4.33) / \Dv-iDv)Qjdxdt<c(^Y -f \Dv - {Dv)q,\ dx 



dt 



holds whenever Qg Q Qr is a parabolic cylinder with the same vertex of Qu. More 
in general, the inequality 



(4.34) 



/ \D^v-iD^v)Qjdxdt<c(^Y J \D^v-{D^v)Q,\dxdt 
Jq, JQr 



holds whenever ^ € {1, . . . ,n} and g < R. 

Proof. We recali the difFerence quotients argument developed for instance in [ITl 
Chapter 8]. We indeed have that v is higher difTerentiable in a smaller cylinder: 
V e L^{to - i?^/4, to; W^''^{B{xo, -R/2))) and moreover, whenever ^ e {1, . . . , n} we 
have D^v € C°(to - io; L'^{B{xo, R/2))). Finally, the function D^v solvcs the 

following diffcrcntiatcd cquation: 

{D^v)t - div {à{x, t)DD^v) = , 

weakly in (3(a;o, to; ^/2) = (5i^/2, where the matrix ày(a::, t) :— Dz.ai{xo,t, Dv{x,t)) 
has measurable entries. By assumptions (|1.18p we have that the following mono- 
tonicity and growth conditions are satisfied: 

i'\z\'^ < {à{x,t)z, z) , \à{x,t)z\ < L\z\ 

whenever z G K", and {x,t) G Qr/2- Therefore, if we let b{x,t,z) := à{x,t)z, we 
may apply Proposition 14. Il to ù = D^v in the cylinder (3^/2 1 thereby obtaining 

/ \D^v-{D^v)Qjdxdt < c(-^) / \D^v-{D^v)Q^^Jdxdt 

(4.35) < \D^v-{D^v)Qjdxdt 

that holds whenever g < R/2, and where again c = c(n,v,L)] we noticc that the 
last estimate in (j4.35|) has been obtaincd arguing as for (j3.43p . Finally, the same 
inequality foUows for g £ (i?/2,i?], trivially, and therefore (|4.34p is completely 
established. In turn, since ^ G {1, . . . , n} is arbitrary (|4.34p implies (|4.33p and the 
proof is complete. □ 
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4.3. Proof of Theorems [HMHH We start by Thcorcm O In Sections OlO] 
we have built the necessary set up in order to adapt the elhptic proof of Section 
13.41 to the parabohc case, and we shaU therefore outUnc the relevant modifìcations. 
Conibining esthiiates (j4.29p and (|4.33p as done m Lemma [3T6l with cstimates p.6p 
and (|3.37p . we have the analogue of estimate (|3.4ip 

/ \Du^iDu)Qjdxdt<ci(-^Y -f \Du~{Du)Q,^\dxdt 
Jq„ Jq-,„ 



QlR 

(4.36) ^^'^l^lf + ^2 (^^y [LMR)]j^ {\Du\ + s)dxdt, 



g j n"^^ \ g 

where the constants c, ci,C2 > 1 and /3 e (0,1] depend only on n,iy,L; the last 
inequahty holds whenever Qg C are backward parabohc cyhnders sharing the 
same vertex, and with R < Ri where Ri = i?i(Li, w(-)) > being a suitable small 
radius; the rcstriction on R is not necessary whcn a(-) is indepcndent of x, although 
a(-) can be stiU depending on the time variable t. By (|4.36p . choosing H G N - 
similarly to (j3.45p - large enough to have ci/H^ < 1/4, and taking R < R < Ri 
such that C2{2H)^[LiLj{R)] < 1/4 similarly to (piig]) . and proceeding as for f^M\i - 
([Og)) we obtain 

/ \Du-{Du)Q^,„\dxdt<]-j \Du-{Du)Q,Jdxdt 

•^Qb/H J Q-2R 

(4.37) [\{Du)q,,\ + s] . 

The last inequality is again valid for cylinders with same vertex and for a Constant 
c depending on n, v, L, with the restriction R < R. We now just have to foUow the 
scheme after p.49p . Specifically, we let 

Q, ■.^Q{xo,to;R/{2Hy), h := \{Du)q,\; 

letting - as in (j3.52p - 

Ai-.^-l \Du~ {Du)Q^\dxdt 

we apply (j4.37p with Qji = Qi, and itcrating as after (j3.53p . we finally achieve 
estimate (|1.23p . 

As for Thcorcm 11.31 to obtain p.25p we argue as foUows. We start by (|4.34p : 
this holds directly for w defined in (|4.8p - (|4.9p . recali that in this case w = w as we 
do not have x-dependence and we do not have to freeze at xq. Next we combine 
(|4.34p directly with (j4.10p and obtain the foUowing analogue of (|4.36p : 

/ \D^u- {D^u)Q^\dxdt < Ci (^) -f \D^u- {D^u)Q^^\dxdt 

JQ, JQ2R 

(4.38) ^JRVWQ^.) 



Now (|1.25p foUows by iterating as in the proof of Thcorcm 11. 21 by setting this time 
g, :^Q{xo,to;R/{2Hy), h := \{D^u)q^\ 

and 

At ■■= f \D^u - {D^u)Q.\dxdt . 
Worth remarking diffcrcnccs are that in (|3.60p the valuc of it is now 
M:^f {\D^u\+ s)dxdt + l'^{xo,to;2R) 

JQixo.to-.R) 
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and this comes from the fact that, being now uj{-) we also have d{2R) ~ and 
therefore (|3.6ip can be replaced by 

Ao + ko<c-l {\D^u\+ s)dxdt. 

jQ{xo,to;R) 

4.4. Proof of Theorem 11.41 The main difference with respect to the proof of 
Theorem ll.2l is that we just need to use estimate (j4.1ip . and no othcr intermediate 
comparison inequaUty of the type (j4.29|l . This is essentially the reason why Theorem 
11.41 apphcs - as expected for zero-order estimates - to operators with measurable 
coefficients. As for the proof, we first obtain the estimate 



/ \u- {u)Q^\dxdt < Ci (-|) / \u-{u) 



2R 

2R. 



dx dt 



(4.39) +c^-j ^^^+ci?. 

which is the analogue of (|4.38p and of (|4.37p when no space coefficients come into 
the play. As usuai, (|4.39p works for backward cylinders with the same vertex, and 
the constants c, ci > 1 and (ì G (0,1] depending only on n,v,L. Estimate (|4.39p 
can be derived by first applying Proposition l4. Il with ù = w and ò(-) = a(-), thereby 
getting the reference estimate 

(4.40) / \w~{w)Qjdxdt<c(-^Y^'^ [ \w~{w)Q^^\dxdt + cRs, 

JQe •'Q2R 

whenever Qg C Q2R is a cylinder with the vertex of Q2R, and where c > 1 and 
P € (0, 1] depend only on n, ly, L. Combining this last estimate with (|4.1ip in the 
same way of Lemma 13.61 - and exactly for (|4.36p - finally yields (|4.39p . By (|4.39p , 
choosing H to have ci/H^ < 1/2, we obtain, again for c = c{n, ly, L) 

(4.41) / \u-{u)Q^,^\dxdt<\i \u-{u)Q,,\dxdt+ '''^lf^^K cRs. 

JQr/H ^ JQ2R ^ 

The last estimate can be iterated as done for (j4.37p . as after (|3.49p : we let 

Q^ :=Q{xo,to;R/{2Hy), h := |(u)q.| 

and finally 



Jq, 



A: := 

and apply (|4.4ip with Qr = Qi, iterating as after (|3.53p . thereby concluding the 
proof. 

Remark 4.1 (Zero order elliptic estimate). The method of proof used for Theorem 
11.41 easily leads to an alternative proof of the zero order estimate (|1.14p - even with 
7 = 1- from [351 ES] for the elliptic case (|l.ip under the assumptions 

\a{x,z)\<L{\z\ + s)P-^ 

v\z2 - Zi 1^ < (a(a;, Z2) - a(x, zi), Z2 - zi) 

whenever zi, Z2 € K" and a; G iì; compare with (|1.28p . We observe that this proof 
works also in the case of general signed measures, and the final outcome is 



(4.42) 



(4.43) \u{xQ)\<cj (|w| +i?s) da; + cW5',p(a;o,2i?). 

Jb{x(,,R) 

The proof in question seems to be slightly easier that those proposed up to now, 
although the assumptions (|4.42p are slightly stronger than those in [SS]. For the 
proof, foUowing the scheme of the proof of Theorem 11. 41 we essentially need to have 
two ingredients: a reference decay estimate for comparison function w introduced 
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in (|3.23p . and a suitable comparison estimate between u and w. This last onc has 
already been proved in p.25p . while foUowing for instance [IHl Theorem 5.1], and 
using as a starting point the proof of the decay estimate in [24l Theorem 7.7], we 
have that the foUowing estimate holds whenever Bg C i?2i? is a concentric to B2R 
for the function w in (|3.23p : 

/ \w - {w)BAdx <c(^ -f \w - {w)B2n\dx + cRs . 

At this stage combining the previous estimate with (|3.25p in the same way as we 
combined (j4.40p with (|4.1ip in the proof of Theorem 11.41 and iterating as in the 
proof of Thcorcm ll.il we finaUy come up with (|4.43p . We just mention that, with 
the notation introduced in the proof of Theorem 11.11 here we have to define 



A;j:=|(m)bJ, and := f \u~{u)B.\dx, 

J Bi 

while no assumption other than measurabiUty is to be assumed on the partial map 
X f-> a{x, ■). 



5. General weak, and very weak, solutions 

5.1. General elliptic problems. Here we consider Dirichlet problems of the type 
(|1.7p assuming that /i e M{VÌ) is a general Radon measure with finite total mass, 
and prove that estimate (|1.6p holds for general weak and very weak solutions. We 
are hereby considering the case of zero boundary datum for simplicity, but more 
general cases can be considered; see also Theorem 15.21 bclow for general weak (en- 
ergy) solutions. We need some preliminary terminology. As mentioned in Section 
2 uniqueness of general very weak solutions fails, thereforc one is led to consider 
very weak solutions enjoying additional propcrties, which in some situations are 
unique; one of such classes, of interest here, is the one of Solutions Obtained by 
Limit of Approximations (SOLA). For the sake of completeness we here recali 
the approximation procedure in question; the main reference here are the works of 
Boccardo & Gallòuet [TUl E] and Dall'Aglio [T3]. We consider a standard, sym- 
metric and non-negative moUifier € C^{Bi) such that ||(/>||li(R") = Ij and then 
define, for every positive integer /i, the moUifier (t)h{x) := h"(l){hx). FinaUy the 
functions iih '■ R" ~^ K are defined via convolution, fih{x) := (/i * (l>h){x). Next, by 
standard monotonicity methods, we find a unique solution G WQ'^{fl) to 

, , r -div a{x, Duh) ^ fJ-h in ri 

^ ' ^ I Uh =0 on do,. 

Up to passing to a not relabeled subsequence we may assume that fi^ ^ fJ- weakly 
in the sense of measures, while the results in (THl E] imply 

71(7) — 1 ) 

(5.2) Duh Du strongly in L'alvi) for every q < , and a. e. 

n — 1 

so that (|1.7p is solvcd by u in the usuai distributional sense, and therefore u is a 
SOLA to (jl.7p ; this of being a limit solution of more regular solutions indeed defines 
SOLA. Moreover, by [Il [15], in the case ^ € L^iyì) we also have that u is the only 
SOLA of (|1.7p . in the sense that if w e W^o'P"^(17) is a distributional solution to 
(ll.7P i obtainable as a pointwise limit of solutions Vh € Wq'^{Vl) to problems of the 
type (|5.ip with fjLh replaced by jìh and jlh ^ weakly in i^(ri), then we have 
u = v; for more Information on uniqueness see Remark 15.11 below. Now we turn 
back to the proof of Theorem 11.11 and in particular to estimate p.55p , that we 
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write when applied to Uh, that is 



(5.3) < c 



A 



ih) 



^0 



E 

1=0 



-c^[Lic.(iì,)]t(fcf)+s), 



i=0 



where here we have obviously set 



\{Duh)i 



and 



A 



ih) ._ 



\Duh - {Duh)Bn]dx 



compare with p.50p and (|3.52p . respectively. We note that the apphcation of 
estimate (|3.55p to Uh is legai sinee the only thing we used in the proof of (|3.55p 
was that u G Wi^^{fl); anyway in the present setting we have Du E C°(ri) by 



standard regularity theory. Letting h 
convergence of measures, yields 



(5.4) k„,+i <clAo + ko+Y^ 



and using (|5.2p and the weak 

m— 1 

+ cY.'^LiLo{R,)]Ì{h + s), 

i=0 



where we are adopting the notation established in the proof of Theorem ll.il From 
(|5.4p on the rest of the proof foUows as for Theorem 11.11 after (|3.55p , but taking 
into aceount essentially two facts: first, the eonvergenee in (j3.66p just takes place 
at Lebesgue points, and therefore almost everywhere; seeond a slight - but obvi- 
ous - adjustment (enlargement of the balls) has to be made in order to overcome 
the presence of \^\{Bi) rather than |/i|(i?i), and to reeover the definition of Wolff 
potential. We have therefore proved the foUowing: 

Theorem 5.1. Let u e W^'^^"^ {Vi) be a SOLA to the problem (|1.7p - which is 
unique in the case /i e L^{Vl) - under the assumptions p.2p and (jl.3p . Then 
there exists a non-negative Constant c = c{n,p,v, L), and a positive radius R = 
Rin^pjh', L, Li,uj{-)) such that estimate (|1.6p holds whenever B{xq,2R) C n and 
R < R, for almost every xq S fi. Moreover, when the vector field a(-) is independent 
of X, estimate (jl.6p holds without any restriction on R. 

Remark 5.1. As proved in [71 [57], in the case p = 2 the Solutions Obtained by 
Limit of Approximations (SOLA) of previous theorem are unique in the sense that 
is a distributional solution to (|1.7P i obtainable as a pointwise 
limit of Solutions Vh S Wq''^~^{ìì) to problems of the type (|5.ip with replaeed 
by fifi and fih then we have u = v. Moreover, as remarked a few lines above, 

the same uniqueness result holds when ^ e L^(ri). 

In the case /i G VK"^'^ (fi), we can deal with standard weak solutions, i.e. u S 
^loc' (^) ' niean of locai approximations of the type 

-div a{x, Duh) = Uh in fi' 

Uh ~ u on dfl' 



(5.5) 



whenever f2' (g fi are smooth sub-domains, just along the lines of the proof above. In 
this case the uniqueness of solutions obtainable in the limit is obviously guaranteed 
as ali solutions are of class W^'P{ìì'). We recali that a characterization of those 
measures fi such that /.j G W~'^'P (fi) is that 



(5.6) 



Wl^ix,l)d\,,\ix) 



n Jo 



\fi\iBix,Q)) 



— d\ii\{x) < oo 
Q 



as proved in [55|; see also [SOI Theorem 4.7.5]. In particular a measure /x G M.{^) 
satisfying the density condition 
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for some e > satisfies (|5.6p : for this case see also [3U[531lin]- We again recali that 
here we have trivially extended fi to the whole K". Moreover the approximation 
scheme in ()5.5p is not necessary when n is an integrable function. We recali that 
the condition /i £ _L"p-"+p (fi) ioi p < n implics fi S W~^'P (il) simply by Sobolev 
embedding theorcm. In this case the proof of Thcorcm 11.11 in Scction [3] works 
dircctly and Icads to (|1.6p which thcrcby holds almost everywhere. Summarizing 

Theorem 5.2. Let u € he a locai weak solution to (jl.ip . under the as- 

sumptions (jl.2p and (jl.3p . and with fi g W~^'^ (fi). Then there exists a Constant 
c = cinjp^v, L) > 0, and a positive radius R = Riji^p,!^, L, Li,lu{-)) such that es- 
timate (jl.6p holds whenever B{xo,'2R) C Q and R < R, for almost every xq G Q. 
Moreover, when the vector field a{-) is independent ofx, estimate (|1.6p holds with- 
out any restriction on R. 

Remark 5.2. By recent higher difFerentiability result for solutions to measure data 
problems - see [481 Theorem 1.5] - in the last two theorems estimate (|1.6p does not 
only hold for a. e. Xq € £1, but indccd outside a set with Hausdorff measure Icss 
than n — 1. In the case fi € Lj^^(ri), with 1 < 7 < np/{np — n + p), the Hausdorff 
dimcnsion decreases up to n — 7. 

5.2. General parabolìc problems. Here, along the lincs of the previous scction, 
we give the extcnsion of estimatcs (|1.23p to parabolic SOLA of p.26p for general 
measures fi € M.{^t)i foUowing this timo the approximation schcmcs in pillSj. To 
this aim we shall consider Cauchy-Dirichlct problems of the type 

I dfUh - div a{x, t, Duh) = fih in 

solved by a functions u,, G C"{~T, 0; L^{fl))nL^{-T, 0; W^-^{fl)), heN. As in the 
cUiptic case fifi is a scqucncc of smooth functions obtaincd from fi via moUification, 
and therefore such that fi^ fi weakly in the sense of measures. The main con- 
vergence result in [5] claims that, up to extracting a not-relabeled subsequence, we 
have that there exists a function u G L'^{—T, 0; W^ '^{iì)) such that Uh ^ u weakly 
in L'"(— r, 0; iy^''(r2)), for every choice of the parameters 

Ti ^ TI 

(5.8) q< -, re [1,2], _ + + 

ri — 1 r q 

see the assumptions of |9l Theorem 1.2], while the result is proved in [9l Section 
3]. In particular we can choose q — r > 1 stili matching (|5.8p and therefore 
the main convcrgence result in [21 Theorem 3.3] statcs that, always up to a not- 
relabclcd subsequence, Du^ Du a. e. in ^t- Then, sincc the sequencc {Duh\ 
is stili boundcd in L'^(ÌÌt) we infer that Duh Du strongly in L^{VIt), and 
finally u weakly solves (|1.26p in the sense of (|2.ip : in othcr words w is a SOLA to 
(|1.26p . We recali that uniqueness of SOLA stili holds in the parabolic case provided 
fi £ L^{ì1t)', for this we refer to [T5]. At this stage, we can repeat the strategy 
of the elliptic proof of Theorem 15.11 i.e. we go to the proof of Theorem 11.21 we 
apply it to every solutions Uh in order to have inequalities as in (|5.3p . we pass to 
the limit with respect to h G N the resulting estimates on the numbers km^ = km 
- compare with (|5.3p - and then we conclude with the desired pointwise inequality. 
Summarizing 

Theorem 5.3. Let u E L^{-T,0;W^''^{n)) he a SOLA to the prohlem (fT:26)) - 
which is unique in the case fi E L^{flT) - under the assumptions (|1.18p and (jl.3p . 
Then there exists a Constant c = c(n, v, L) and a radius R = i?(n, v, L, Li, lu{-)) such 
that estimate (jl.23p holds whenever Q(xo,tQ;2R) C fi, for almost every (xqjÌq) G 
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ÌÌt, provided R < R. When the vector field a(-) is independent of the space variable 
X, estimate ()1.23p holds without any restriction on R. 

We similarly have 

Theorem 5.4. Let u e L^{-T,0;W^^\n)) he a SOLA to the problem - 
which is unique in the case fi G L^(Q,t) - under the assumptions (|1.28p . Then 
there exists a Constant c = c{n, L) and a radius R = R{n, L, Li,lu(-)) such that 
estimate (|1.23p holds whenever Q{xo,to;2R) C fi, for almost every {xojÌq) € Qt, 
provided R < R. When the vector field a(-) is independent of the space variable x, 
estimate (jl.23p holds without any restriction on R. 

Morcovcr, when fi is regular cnough - for instancc when it belongs to the dual 
space T, 0; W~^'^{n)) - estimate (|1.23|) works for general weak solutions u G 

6. INTEGRABILITY ESTIMATES 

In this section we rapidly show how to get the integrabihty cstimates directly 
from pointwise cstimates as (jl.6p and (|1.23p ; the most interesting results are prob- 
ably those in the parabohc case, since they are completcly new, as described below, 
while in the eUiptic case we shall mostly find a imified approach to the proof of 
several results appearing in the literature. In both the elliptic and the parabolic 
case we shall confine ourselves to outline the strategy of proof, since then the details 
can be easily added by the interested reader. 

A main ingredient in the elliptic case here is of course the possibility to control 
the non-linear Wolff type potentials via the non-linear Havin-Maz'ja potential, and 
ultimately via Riesz potentials. which is described in (jl.lSp . Therefore the main 
property to use is the boundedness of the Riesz potential in Lebesgue spaees 

(6.1) <c||5||l.(B") ge(l,7i//3), 

with similar cstimates to hold also in Marcinkiewicz and Orlicz spaees via inter- 
polation; see also |35] and next section. Using (|6.ip and (|1.15p we then come to 
various mapping properties as for instancc 

(6.2) fieL^' =^W'i i;R)eL^'^ ge(i,n). 

p 'P 

The last implication comes of course with related quantitative cstimates. In turn, 
using (|6.2p in combination with (|1.6p leads to deduce the foUowing inclusion for 
solutions to (|1.7p considered in Theorems l5.Hl5T2l 

(6.3) fi e =^ Du e (fi) qe{l,n), 

The last implication allows to recover various cstimates available in the literature 
pn rm irsi which are usually achieved via very different techniques, according 
to the size of q. There is now a large number of possible variants to the last result. 
In fact by the known bounds on the Riesz potential, it is possible to easily find 
virtually ali kinds of rearrangement invariant function spaees integrai cstimates. 
For more details on this we refer to the next section, which is dedicated to parabolic 
problems, and to Remark l6.1l below. Further developments of our viewpoint can be 
found in the forthcoming paper of Ciancili [Mj , where a few interesting consequences 
of estimate (|1.6p are presented. For instance, in [T3] the following consequence of 
(|1.8p is presented 

(6.4) fiE L[n, -^—] =^ Due locally in fi . 
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In this case something more can be actually asserted, namely the continuity of 
Du, as will be proved in the subsequent paper [53]. Again, Wolff potentials based 
continuity criteria for the gradient of solutions can be obtained starting from the 
methods presented here, and will be presented in [23j . 

We also mention that by localizing inclusions (|6.2p and (|6.3p it is also possible 
to obtain explicit locai cstimates rclated to the rcgularity results mentioned above, 
something which is also not casily reachable via the global techniques usually used 
in the literaturc. For this we also rcfcr to P71H51H5] . 

Remark 6.1 (Borderline cases). In the regularity theory of measure data problems, 
for a long while it has been an open issue to establish the validity of estimates in 
Marcinkiewicz spaces in the conformai case p = n. Specifically, for a SOLA solution 
u of (|1.1|) we bave the inclusion Du S Al" (fi) - see (j6.10[) below for the definitions. 
This has been first proved in [IH] directly for the p-Laplacean system, see also [17] 
for different approaches. Here we would like to outline how estimates (jl.6p - (|1.16p 
immediately imply that Du € AljJ,^(0,M") in our sctting. This goes via analyzing 
properties of Wolff potentials in Lorentz and Marcinkiewicz spaces - see (j6.9p - (j6.10p 
below for the definitions; indeed by standard properties of Ricsz potentials - see for 
instance [IHl Section 4] - and (jl.lSp we bave 

More in general, as in (|6.2p we have that (|1.15p iniplies 

(6.5) M e L{q,j) =^ wl (,i?) e L ('!^fcll,^(p- 1)) 

whenever q € (1,?^) and 7 G (0,cx)]. Anothcr open regularity problem was the 
borderline L^-regularity of solutions of SOLA to (|1.4p . this means cstablishing, in 
the case p < n, that 

(6.6) Ai e L ( ^\=^Du(^ L{p, j{p - 1)) 

\np — n + p J 

whenever 7 S (0, 00]. This problem, raised several times in the literature - see for 
instance [SI [31] and related references - has been recently settled in [JS] Theorem 
2]. Here another proof of (j6.6p foUows directly from estimates (|1.6p - (|1.16p via (|6.5p . 

Remark 6.2 (Sharpness). Inequality p.6p is in a scuse optimal as W'^ allows 

to recast ali the integrai estimates in the scale of Lebesguc and Lorentz spaces as 
describcd in the prcvious remark. This is obvious in the case p ~ 2 due to well- 
known representation formulas for the Poisson equation (|1.12p : the optimality for 
p > 2 foUows since the result in (|6.3p is the best possible; at this point a better 
potential in (|1.6p would imply a better estimate in (|6.3p . which is indeed impossible. 
In other words, in (|1.6p no potential of the type with (3 > 1/p can replace 

, for every p > 2. Lct us now test estimate (jl.6p in the case we have the 

p 'P 

following Dirichlet problem for p < n: 

Sn in Bi 
on dBi 

with Sn := CT„_i(5, where S is the Dirac measure charging the origin and ct„_i := 
?i"~^(9i?i) is the (71 — l)-dimensional Hausdorff measure of dBi. The unique 
solution to (|6.7p is given by the so called fundamental solution of the p-Laplacean 



(6.7) 



—ApU 
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operator, that is 

(^\x\^ -ij if p <n 
^ — log |a;| ii p = n . 



u{x) 



The uniqueness is a consequeiice of classical results of Serrili [SI] and Kichenassamy 
& Veron [3T]; an account of the proof can be found in [SOI Section 4.4]. Let us show 
that in this case estimate (|1.6p reverses: let us take xq such that |a;o| < 1/2, 
then we have 

(6.8) / \Du\dx + W{" {xo,2R)<c\Du{xo)\, 

with R ~ \xq\, and for a Constant c depending only on n and p. Indeed by co-area 
formula we have 



f \Du\ dx < c{n) 4 

2R p 

\Du{y)\dH"-\y)dQ 



\Du\ dx 

lB{xaM) J B{0,2R) 

I \ r2R 
cin) 



c{n,p) f'^^' „+i^ dg 

= — / Q — 

ii" io Q 

= c{n,p)R~p^ < c{n,p)\Du{xo)\ . 

On the othcr hand, as 5n{B{xo, g)) = for g < R and 5n{B{xo, g)) = 1 othcrwisc, 
we also have 

W*" (xo,2i?) < c{n,p) / g^ — < c{n,p)RT=r < c{n,p)\Du{xo)\ , 
v'P Jr g 

so that (|6.8p foUows combining the last two incqualities. Inequalities as (|6.8p do not 
only hold when the right hand side of the equation is a Dirac measurc; indeed, when 
considcring more general measures concentrating on lower dimensionai set - lines, 
hypcr-plancs, etc - we observe powcr-type singularitics of solutions and cstimates 
like dea. 



6.1. Parabolic integrai estimates. In the parabolic case the consequences of 
the pointwise estimates are more interesting as we are able to get sharp borderline 
estimates in Lorentz and Marcinkiewicz spaces which do not seem immediately 
reachable via the known parabolic techniques. Moreover, as in the elliptic case, 
we catch directly borderline cases in Lorentz spaces. To this aim let us recali that 
a map g : A C R^ ^ with A being an open subset and fc S N, and for 

q G [1, oo) and 7 e (0, 00], is said to belong to the Lorentz space L{q, ^)[A) iff: 

(6-9) ll.9lll(,,,)(^):=7y^ (A'|{.TeA : |g(x)|>A}|)? -<cx) for 7 < ex. . 

In the case 7 = 00 we have instead Marcinkiewicz spaces Ai'^{A) = L{q,oo){A): 
(6.10) sup A«|{x G A : \g{x)\ > X}\ =: ||.9||^,(^) < 00 . 

A>0 ^ ' 

The locai variants of such spaces are then defined in the obvious way. We bere just 
recali that Lorentz spaces "tune" Lebesgue spaces via the second index according 
to the foUowing chain of inclusions: 

= L{r, r) C L{q, 7) C L{q, q) C L{q, r) C ^(7, l) = L\ 
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which hold whenever 0<7<(7<r<oo. The foUowing inequalities can now be 
obtained using the definition of calorie Riesz potentials given in (|1.20p : 

(6-11) ll-^/3(.9)|li(_j^,^)(R.+i) < c||5|U(5,^)(R.+i) 

which holds whenever q > l and f3q < N, and, in the case q ~ l, 

(6.12) ||/^(g)|| « <c\g\iW"+^). 

The last inequaUty stili holds when g is a measure and (3 < N . Inequalities (|6.1ip 
and (|6.12p can be localized, more precisely it is not difhcult to see that //3(-) controls 
the localized potential R) defined in (jl.22p in the sense that R) < Ipig), 
holds. Therefore we have that 

(6.13) I^(.;i?):i(g,7)^i(^_^,7^ and I^(-; i?) : ^ Al^ 

hold, continuously, assuming (3q < N and (3 < N , respectively. The estimates in 
(|6.13p refer to M"+^ as supporting domain, but they can be localized using cut-off 
functions and therefore yield locai estimates for calorie Riesz potentials; see for 
instance [ISlUn]- The final outeomc is the foUowing result: 

Theorem 6.1. Let u e L^{-T,0;W^^^{n)) be a SOLA to the problem (fL^ - 
which is unique in the case ^ G L^{flT) - under the assumptions p.lSp and (jl.3p . 
Then the foUowing implications hold: 

(6.14) ^ e L{q, 7) =^ Du e L ^ ^'^ ,7^ for 1 < q < N and < 7 < 00 
and 

(6.15) e M{n) =^ Du e 

with ali the previous inclusions being meant as locai in Q. 

The previous result is of course a straightforward consequence of estimate (jl.23p 
together with the properties in (|6.13p . Moreover, using a localized version of the 
estimates relative to (|6.13l) it is also possible to obtain explcit locai estimates relative 
to the implications in (j6.14p - (j6.15p . See for instance the strategies adopted in [IS]. 
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